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I.  GLOSSARY  OF  SYMBOLS  USED 


P..  =  correlation  coefficient  between  variables  X.  and  X. 
U  1  J 

PrO  =  probability  of  an  event  occurring 
I  =  correlation  matrix 

X  =  n  dimensional  random  variable  'vN(0,l) 

Xi  *  i-th  variable  component  of  X 
=  Pr{X^  (-z  ,z) } 


A  =  Pr  {  u  (A. ) } 
i-1  1 

a  =  Fixed  upper  bound  for  Pr{A} 

F 

0^  =  Upper  bound  for  Pr{A}  for  a  given  zp 

zF  =  Fixed  edge  length  of  cube  for  =  Pr {X^ el  (-z^. ,  Zp) } 

zQ  =  Edge  length  of  cube  for  *  PHX^C-Zq.Zq)  }  determined  so  that 
Pr{A}  £  aF 

T  =  Tree  connecting  n  disjoint  events 

Tq  =  Best  bivariate  tree  connecting  n  disjoint  events  such  that 
n 

. E1Pr(A.)  -  . Pr(A.nA.)  is  minimized 
1=1  i'  ijeTq  i  3 

B  .  =  {A.uA.}  where  A.  is  linked  to  A.  on  Tn 

13  i  3  1  J  0 

T’  =  A  tree  linking  B..  formed  by  linking  B..-B..,  iff  i=i’  or  i= j ’ , 

7  ij  13 

or  j  =  i ’ ,  or  j=j  ’ 

<p  -  Standard  univariate  normal  density 


=  Bivariate  normal 


density  of  j]] 


ff0l  f1  "l  »2l] 


?3P1’P2,P3  =  Trivariate  normal  density  of  N  10  ,  p^  i  p^ 

•  0  ^p2  p3  1 


iVii 
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DERIVING  AND  APPLYING  IMPROVED  UPPER  BOUNDS  FOR 
MULTIVARIATE  NORMAL  PROBABILITY  OUTSIDE  OF  N-CUBES 


II.  Introduction 

The  probability  areas  outside  n-dimensional  multivariate  normal 

cubes  have  long  been  of  interest  to  statisticians.  If  X  is  a  random 

vector  with  an  n-dimensional  multivariate  normal  distribution  with 

u.  the  mean  of  X.  and  c.  the  standard  deviation  of  x., 
xxi  1 

i  =  1, . . . ,n,  then 


(1)  Pr  X  is  not  in  the  multidimensional  rectangle 


i.-z-z  <  X1  <  u -+z_o-  ' 
1  F  1  1  1  F  1 


V . — 2_0 .  <  X.  <  U.+2_0. 
i  F  i  1  1  F  l 


U  -ZrG  <  X  <  u  +z_o 
-  n  F  n  n  n  F  n 


is  the  same  as 


(2)  Pr  Y  is  not  the  multidimensional  cube 


X.-y. 

where  Y.  =  — — — 
l  a . 

l 


i.e.  Y  is  X  witn  each  component  standardized  to  N(0,1). 


rzF  <  yi  <  zf1 
-2f  <  <  zF 

-ZF  "  Yn  "  ZF 


This  probability  is  of  interest  in  three  cases: 

A.  Finding  the  probability,  that  a  random  observation  X  is 

not  in  the  region  determined  by  (1)  when  u  ,u_,...,u  are  known. 

1  2  n 

B.  Simultaneous  confidence  interval  construction  producing  confidence 

intervals  for  u.  from  X.  of  the  form  [X.-z  o.<p.  <  X  +z  n  1 
i  i  l  0  l  l  l  0  l 

and  being  at  least  (l-a_)  confident  that  all  confidence 

intervals  cover  the  true  parameters.  This  is  done  when  u  ,u  , . .  .  ,u 

12  n 

are  unknown. 

C.  When  simultaneously  testing  the  hypothesis  H^:  * u  , 

H,,:  :U  =  u  against  the  alternative  hypothesis 

22Znnn 

H,  :H..  is  false •••H:  H  is  false,  and  it  is  desired  to  obtain  an 
j_a  ±  na  n 

upper  bound  otp  for  the  probability  that  at  least  one  is  rejected 
when  in  fact  all  are  true  by  using  acceptance  regions  constructed 

the  same  way  as  in  case  B. 

In  case  A,  is  usually  fixed  in  advance  and  it  is  desired  to 

determine  the  smallest  upper  bound  for  a.  In  cases  B  and  C,  Up 

is  usually  fixed  in  advance  and  it  is  desired  to  determine  the  smallest 
Zq  which  will  give  Op  as  an  upper  bound  for  the  true  a. 

Current  methods  used  to  find  and  z^  in  the  situations 

described  above  are  of  three  types  (i)  the  Scheffe  method  which  incor¬ 
porates  correlations  of  marginal  distributions  along  the  axis  of  the 
cubes  but  builds  the  confidence  intervals  for  all  linear  combinations  of 
the  marginal  axis  variables,  and  this  gives  extremely  conservative 
(large)  values,  (ii)  the  Tukey  range  and  Bonferroni  inequality  methods 
which  generally  give  less  conservative  results  but  do  not  incorporate 


1 

a 


| 


>  A  .V.VAV-'V.V.'-A.'-.V.V 


correlations  of  marginal  distributions  along  the  axis  of  the  cubes, 
(iii)  A  method  suggested  by  Hahn  uses  the  (not  necessarily  correct) 


assumption  that  if  P _  _>  0  for  all  i,j  then 


...  *ntt)dx  ...dx  > 

-z  1  -z  1 -z 


$  (P)dx  . . . dx  , 
n  1  n 


where  (a)  Z  is  the  covariance  matrix  of  x,  ,...,x  .  WLOG  let  the  variance 

1  n 

*  . 

of  all  x.  be  1,  the  diagonal  elements  of  Z  are  1.  (b)  P  =  P 

■k  |  | 

where  P  is  the  minimum  over  i  j  of  iP..  the  absolute 

1  ij  1 

correlation  between  x.  and  x..  (c)  $  (P)  is  an  n  dimensional 

i  3  n 

multivariate  normal  distribution  with  var(X.)  =  1  for  all  i  and 

l 

cov(x.x.)  =  P  for  all  i,j.  A  problem  with  this  method  is  that 
often  when  n  is  large,  P*  is  likely  to  be  very  close  to  zero  giving 
this  method  little  advantage  over  a  method  assuming  independence. 

The  method  described  in  this  report  is  an  improved  version  of  the 
Bonferoni  method  which  can  incorporate  correlations  of  the  marginal 
distributions  to  produce  even  less  conservative  estimates  for  Zq  and 
Uq  than  does  the  standard  Bonferroni  method,  or  any  of  the  other  methods 
mentioned  above. 

Ill .  Description  of  the  Improved  Bonferroni  Method 

a.  Three  Pathway  Overlap  Reduction  Theorem.  Hunter  (1976),  Worsley 
(1982).  Represent  events  A^,...,An  as  vertices  V^,...,V  of  a  tree 
whose  vertices  and  V  are  connected  are  joined  by  edges  e^ , 


Pr  U  A  <  l  Pr (A. ) 
i=l  i=l 


L  Pr(A  r,A  ) 

i,j:e..cT  J 

ij 


|u 


I 


1 

Eh 


#5(*i 

W! 


W3 

w!l 

i 


E 


Proof :  Since  T  is  a  tree,  it  is  always  possible  to  find  a  permu¬ 


tation  P, of  l,...,n  so  that  A  is  joined  to  A  for 

in  p  p . 

i  J 

some  j  <  i  (i=2,...,n).  We  can  then  write 


(4) 


11 

(Ja  =A  U  (A  ^  A  n  A  ) 
i=l  Pi  Pi  Po  Po  P 


'J  •  •  •  u 


2  F2  f1 


(A  ^  A  n(A  u* • • u  A  )). 
Pn  Pn  Pn-1  P1 


The  events  in  (4)  are  disjoint,  so 


(5)  Pr 


n 

,  U  A 

Li-1  pl 


Pr(A  )+Pr(A  -vA  oA  )+•••  + 
P1  P2  p2  P1 


Pr(A  'V  (A  n  (a  u  ...  uA  )) 
Pn  pn  pn-l  P1 


but  Pr(A  \A  n  (A  j  •••  uA  ))  <  Pr(A  %  A  n  A  )  where  1  <  j  <  i- 
Pi  Pi  Pi  Pi-1  pi  pi  Pj 


(5)  <  7  Pr (A  )  -  1  Pr  (A  "  A  ). 

pi  i,j:e..cT  pi  pj 
ij  J 


Venn  Diagram  Explaining  Theorem 


\  A2  A3  \  A5  A6 


[Note  (A^  and  A^)  (A2  and  A^) ,  etc.  are  not  necessarily  disjoint 


1  hence 


but  are  drawn  that  may  be  due  to  inability  to  represent  6  dimensional 
intersections  on  a  2  dimensional  diagram. J 


This  is  a  tree  since  there  is  no  loop,  i.e. 


A  -*■  B 

t  + 
D  -  C 


6 

So  Pr(  u  A.)  <  Pr(A..  )+Pr  (A„)+Pr  (A^)+Pr  (A.  )+Pr(A^)+Pr  (A,)-Pr(A1  iA.)- 
1  —  i  A  j  4  0  o  i  Z 

Pr(A^  n A^)  -  Pr(A^  n A^)-Pr(A^  n A^)-Pr(A^  n  A^) . 

A  theorem  will  be  presented  in  a  future  technical  report  of  this 

series  written  by  the  same  author  (Hoover  1986)  which  will  enable  us,  in 

n 

some  cases,  to  get  a  lower  upper  bound  for  P  (  A.)  than  did  the  previous 

r  i=l  1 

theorem  from  incorporation  of  probabilities  of  pairwise  intersections  of 

events.  The  improvement  (when  there  is  one)  will  often  be  very  miniscule, 

and  in  practice  this  theorem  will  be  computationally  much  more  difficult 

to  apply  than  is  the  three  pathway  overlap  reduction  theorem. 

IV .  Application  of  Improved  Bonferoni  Method  to  Multivariate  Normal 
Distributions. 

a.  Bivariate  Method  -  Incorporating  Only  Bivariate  and  Univariate 
Normal  Distributions . 

If  [Y]  is  an  n-dimensional  normalized  random  vector  with  mean 
n 

[0]  and  covariance 
n 


1  P 


'12 


12... 


In 


L  ”ln  •" 


1  *. 
. . .  1 


and  z  >  0.  Then  let  the  event 


[V2'V2 


(6.5) 


and  Che  event 


A  =  < 


Y  -z,Y  +z 
n  n 


does  not  include 


(6.6)  =  (Y^-z,Y^+z)  does  not  indlude  [0] 


fy . -z ,v . +Z) 

LEMMA  (I)  (Sidak,  1968)  Pr [ A .  n A J  =  Pr  1  1  '  does  not 

l  j  [V ,-z ,v ,+z i 

0  “  J  '  J  , 

include  [^]  is  a  monotonically  increasing  function  of  j j  I  :  that 


lptjl  >  lptt!  ^PrlAj^A,]  >  Pr[AkTA. 


We  should  also  note  that 


(7.5)  Pr(A.)  =  1  -  j  f(x)dx  for  all  i  =  l,...,n 

2  -z 


where  ■>  is  the  standard  normal  density  and  f  b(x)dx  can  be  computed 


using  the  IMSL  subroutine  MDNOR;  also 


(7.5’) 


Pr  (A,  A  . )  =  Pr  (A . )  +  Pr  (A  .  )  -  PrCA,  j  A.) 
1  J  i  3  1  J 


where  Pr(Ai)  and  Pr(Aj)  are  described  in  (7.5)  above  and  Pr(A^  J A  ) 
can  be  found  using  the  bivariate  normal  distribution 


f  z  z 

(7.6)  Pr(A.  A.)  =  1  -  I  : ■^P\x1,x?)dx1dx0 

3  ' -z  ' -z 


m 

§ 

ft 


can  be  constructed  such  that  if 


over  the  i  and  j 


Theorem.  A  best  bivariate  tree 
N_  =  a  vector  of  the  order  statistics  of  |P,. 

o  ij 

linked  on  then  for  all  k  =  (l,...,n)  and  all  other  trees  T 

it  will  be  true  that  N  (k)  >  N„(k).  This  best  bivariate  tree  T_ 

T0  “  T  0 

is  constructed  by  the  following  procedure: 


(8)  Perform  the  following  (n-1)  times.  Link  together  any  two 


variables  x.  and  x.  such  that 
1  J 


P^.  |  is  a  maximum  among 


all  links  which  have  not  been  established  in  a  previous  step 
and  whose  establishment  will  not  form  a  loop  with  the  other 
links . 


Proof . 

FIRST,  Given  n  points, at  most  (n-1)  links  can  be  made  between 
them  without  forming  a  circuit,  and  it  is  possible  to  have  (n-1) 
links  in  a  tree  no  matter  how  the  tree  started. 

Proof .  Induction 

(i)  Let  n  =  2.  Obviously  only  one  link  is  possible. 

(ii)  Now  if  a  tree  has  been  established  for  m  items  and  another 
item  is  added.  This  item  can  be  added  anywhere  in  the  tree 
either  by  directly  connecting  to  another  element  of  the 
tree  thus  adding  one  link,  or  by  connecting  the  new  ele¬ 
ment  to  P  elements  on  the  tree  and  breaking  exactlv  one  link 
in  each  of  tiie  (P-1)  segments  which  connect  those  P 
elements . 


A  -  D 
I 

B  -  0 


add 


x 


A 

X 

X 

B 


x 

n 


a  -  n 

i 

or  X  or 

B  -  C  X 


A  -  P 
I 

B  -  t: 


8 


km 


The  (P-1)  segments  exist  since  the  P  elements  are  connected 


but  it  is  not  P  or  more  segments  since  there  are  no  circuits. 
Thus  the  net  number  of  new  links  added  to  the  tree  by  bringing 
in  the  new  point  is  P-(P-l)  =  1. 


SECOND 


Now  that  at  most  and  exactly  (n-1)  links  are  possible  in  every 

tree  it  remains  to  be  established  that  N_  (k)  >  N„ (k)  for  all  k. 

T0  ~  T 

This  proof  will  be  inductive.  For  i  =  1  it  is  obvious  that 


N  (1)  >  N„ (1) .  Now  assume  that  for  i  =  l,2,...,k-l  it  has  been 
T0  ~  T 

shown  that  N_  (i)  >N\(i).  Then  assume  that  for  k  there  is  a  tree 


0 


~  T 


T  such  that  N„(k)  >  N  (k) .  Let  J-0  be  the  connection  in  T 
T  Tr> 


which  generates  N„(k).  Then  if  we  consider  the  first,  second,..., 

T 

(k-l)th  links  established  by  the  algorithm  of  this  theorem  in 

forming  TQ  we  must  have  exactly  one  string  connecting  J  and  0, 

i.e.  J  -•••-  L  -  M  -•••-  0.  Otherwise  on  the  kth  step  the  connection 

J-0  is  possible  and  will  be  made  before  a  connection  generating  a 

lower  absolute  correlation  coefficient  is  made,  k-1  links  have  been 

established  by  the  order  statistics.  N  and  k-1  different  links 

0 

(some  of  which  may  be  the  same)  have  been  established  by  the  order 

statistics  N».  The  (k-1)  links  established  bv  the  order  statis- 
T 

tics  N  must  prevent  the  establishment  (excluding  redundancy)  of 
0 

any  of  the  k-1  links  established  by  the  order  statistics  N„.  For 

T 

if  there  exists  a  link  n-C  corresponding  to  the  Pth  order  statistic 
of  N.:  P  <  k-1,  then  |P  1  >  :P  I  >  kth  order  statistic  of 

t  -  X-:’V  ~  xj’xo 

N  and  hence  the  algorithm  on  the  previous  page  will  establish  the 
0 


link  n-5  on  the  kth  step  before  establishing  a  link  with  a  lower 
absolute  correlation. 

LEMMA.  (II):  (k-1)  links  which  do  not  form  a  circuit  prevent  the 

establishment  of  any  of  a  set  of  k-1  different  links  (excluding 
redundancy)  which  themselves  collectively  do  not  form  a  circuit  iff 
each  set  of  the  k-1  links  is  a  collection  of  maximal  spanning  trees 
where  each  maximal  spanning  tree  in  one  set  connects  exactly  the 
same  points  as  one  of  the  maximal  spanning  trees  in  the  other  set. 

Proof  of  Lemma. 

<?=  WLOG  assume  each  set  contains  only  one  maximal  spanning  tree 
and  that  the  maximal  spanning  tree  in  each  set  connects  exactly  the 
same  k  points.  Then  to  add  any  of  the  links  in  L ^  which  are  not 
in  L^  to  the  structure  formed  by  the  links  of  L^  will  create  a 
structure  with  k  links  between  k  points  and  hence  must  contain  a 
circuit.  If  there  are  more  than  one  maximal  spanning  trees  in  L^ 
and  each  maximal  spanning  tree  in  L^  connecting  the  same  points 

as  a  maximal  spanning  tree  in  L^.  The  proof  can  be  extended  bv 
considering  all  the  maximal  spanning  trees  of  L^  and  adding  in  the 
maxima],  spanning  trees  of  one  at  a  time. 

=*>  If  each  tree  in  the  set  of  maximal  spanning  trees  formed 

by  L^  does  not  connect  exactly  to  the  same  points  as  a  corresponding 
tree  in  the  set  of  maximal  spanning  trees  formed  by  L^,  then  find  a 
connection  in  L^  which  connects  two  points  which  are  not  both  in  the 
same  maximal  tree,  formed  by  L^.  Then  adding  this  connection  to  L^ 
will  either  (i)  connect  two  separate  trees  together  which  doesn't  form 


a  circuit  or  (ii)  connect  an  unlinked  point  to  a  tree  which  doesn’t 
form  a  circuit  or  (iii)  connect  two  unlinked  points  together  which 
does  not  form  a  tree.  QED. 

But  this  lemma  creates  a  contradiction  with  the  situation  we 
observe  since  J  cannot  be  connected  to  0  due  to  the  first  k-1  links 
of  Tq  and  each  maximal  spanning  trees  formed  by  the  first  (k-1) 
links  of  Tq  will  contain  exactly  the  same  points  as  a  maximal 
spanning  tree  formed  by  the  first  (k-1)  links  of  f.  Thus  if  it 
were  possible  to  connect  J  and  0  together  on  the  kth  link  of  N- 
then  it  would  also  be  possible  to  connect  J  and  0  together  on  the 
kth  link  of  which  would  have  been  done  by  the  algorithm  forming 

Tn  since  (k)  >  N  (k) . 

0  T  T0 

Implementation  of  the  Procedure  Building  Tq  on  Computer 

This  can  be  done  by  a  program  which  picks  the  maximum  correlation 
coefficient  that  has  not  been  used  for  a  link  and  adds  that  link  to 
the  graph  if  doing  so  does  not  create  a  loop,  i.e.  start  with  F=0  and 
0)  Do  while  F  <  n. 


1)  Pick  largest  remaining  'P^i  i^j . 


P.  .  |  =0. 
1  ij  1 


2)  Link  i  to  i  if  i  and  j  are  not  in  the  same  tree.  Other¬ 
wise  set 

3)  Note  which  tree  i  (or  j)  has  been  added  to  if  a  variable 
has  been  added  to  a  tree.  Otherwise  if  two  trees  have  been 
merged  note  that  the  trees  are  now  the  same. 

4)  F  =  F+l  go  to  step  0. 

3 

This  operation  involves  0(n)  steps  and  for  20  or  less  variables  this 
should  not  be  too  complicated. 


Let  [X]  v  N(0,I),  where  (I)  is  a  correlation  matrix  and  let 


Tg  be  the  best  tree  obtaind  for  X  using  formula  (8) .  If  we 
let  A  be  as  defined  in  (6.5)  and  as  defined  in  (6.6)  for  some 

fixed  >  0 ,  then 


Pr(A)  <_  i  Pr(A.)  -  l  Pr(A.r>A.) 

i-i  1  1  J 


=  nPr (A  )  -  T  Pr(A.  n  a.) 

1  «£To  1  J 


(9)  by  (7.5),  (7.6)  =  n  1  -  $.(x)dx  -  T  j  1  -  $(x)dx 


LC  Zf 


n  r  [zf  fz 

(x)  dx  -  1  -  j 


Zf  rZc  P.  . 

ij1"5  (x1,x2)dx1dx7 


zf  -zf 


=  (2-n)  f  f  i (x)dx  +  V  1 


[Zf  [2f  / 

I  2 

J  _ i?  •  _ 


1^£T0  ~  zf  ~zf 


P  .  . 

]  (x1,x7)dx1dx2 


If  a  is  fixed  in  advance,  say  a  =  a  and  one  wishes  to  determine 

F 


Zq  such  that 


Pr(A|z0)  <_  a?  . 


This  can  be  done  iteratively  and  quite  inexpensively  by  using  the  Newton- 
Raphson  formula.  Start  with  (z 0)  an  educated  guess  for  Zq.  Then  let 


8S 


where  f(Zg)  is  [  (9)  — ctp ]  and 


f,<v  -£i-  ■ 


^i.  (2-n)  1  -  [  $  (x)dx  ■+  l  1-f  f  (x^*x2>dxtcixol 

Zn  *-  J  -Z  -  iiCT-  “  ^  -Z  1-7.  ~~ 


-t\z  1  -  j_  <Kx)dx]  -  -2<Kzq) 


li  i£  \}  '  j  j  i2lj<x1.x2)dxidx]  =  -4Mzq)  f  °  i(x-P  zQ)dx  . 

zo  L  ”z  -z  J  "20 

Therefore 

*  -  ~  rrzn  .  1 

(10)  (*)  =  (2n-4)  <})  (Zq)  -  4?(zq)  )  |jU  $(x-P..zQ)dx  . 

ijETo  J"2o  1J 

Then  set  (Zq)  to  equal  Zq  and  use  the  above  procedure  again  to  produce 
a  better  estimate  for  Zq.  It  turns  out  in  practice  that  the  Newton-Raphson 

A 

method  (z^)  estimate  converges  quite  rapidly  to  z^  and  is  computationally 
very  inexpensive.  Note  that  f  is  monotonic  so  the  solution  is  unique. 


b .  Trivariate  Method  -  Incorporating  Trivariate  and  Bivariate  Normal 


Distributions. 

If  one  is  able  to  integrate  over  trivariate  normal  cubes  (which  can 


L* 
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>K4‘ 
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be  done  by  the  IMSL  subroutine  DMLIN,  but  is  very  expensive  and  time 
consuming]  then  one  can  use  the  method  below  to  obtain  closer  upper 
bounds  for  of  Statements  B,  C  on  page  1  and  of  Statement  A 

on  page  1  than  are  obtained  by  the  Bivariate  Method. 

Let  Tq  be  the  best  bivariate  tree  as  constructed  in  (8).  Then 
for  all  links  ijeT^  define  the  event  B„  =  (A.  uA.),  There  will  be 
(n-1)  events  B  ,  since  there  are  (n-1)  links  on  the  tree.  These 
events  can  then  be  linked  together  in  a  new  tree  T^  where 


(13)  B_  can  be  linked  to  B^,..,  only  if  i=i'  or  i=j  '  or  j  =  i'  or  j=j  1 


For  example  given  events  linked  by  T^  as  below 


\/  A5 


(  A_( 

Ui 

V 

Vv  27/ 

Vv  477 

”6  J 

Then  B12-B23-B24-BA5-B46  is  a  tree,  so  is  B^-B^-B^-B^-B^ . 
Except  in  the  case  where  A^  -  A2  -  •  •  •  -  A^  is  a  string,  there  will 
not  be  a  unique  tree  linking  the  B  '  s  together  such  that  condition 
(13)  is  met. 

Assuming  that  we  have  a  tree  of  events  B  formed  as  in  the 

method  above  with  B„  =  (A^uA^)  where  X  'v  N(0,$)  and  Var(X^)  =  1 

for  all  i  with  A.  =  (|X.|  > z} .  Then 

l  l 


Pr (B . .)  =  Pr (A .  u  A.) 
il  i  3 


ER0 


1 
1 


given  in  formula  (7.6)  and  Pr(B. .  n  B.,  )  =  Pr(B. . )+Pr(B  )-Pr(B. .  j  B  ) 

ij  jk  ij  Jk  ij  jk 

where  Pr(B, .)  and  Pr(B..)  are  as  defined  in  (14)  and 
ij  jk' 


(15) 


Pr  (B .  .  j  B..  )  =  Pr  (A.  uA.uA.) 
ij  jk  i  1  k 

fZ  ( Z  fZ 


=  1  - 


J. 


z  J  -z 


-z 


>3lPij*Pjk»Pjk]dxldx2dsc3 


where 


is  a  3-dimensional  normal  distribution  with 


1 

P .  . 
ij 

Pik 

p. . 

1 

P., 

P1J 

jk 

-  ik 

Pjk 

1  - 

4>-.[P .  .  ,P ,P.,  ]  can  be  claculated  with  DiSL  subroutine  DMLIN  if 
J  ij  ik’  jk 

is  of  full  rank.  Otherwise  there  is  a  procedure  developed  by 
Uusipaika  which  can  be  used. 

After  incorporating  (14)  and  (15)  into  (3),  it  follows  that 


Pr (uA . )  =  Pr(  y 

1  ^  fT' 

■ij  0 


1 


B . .  |  <  7  Pr(B .  .)  -  L  Pr (B . .  n  B., .,) 


U  ”Ti  ij.i  -  ij  4  A  ij  i'j 

•  •  — ■  ^  r\  ■*-  r\ 


0 


0 


L 

To 


1  - 


-z  '  -z 

r. 


)2(Pij)dx1dx2 


ij  *i' 3  ,£TqLL 


-z  '  -z 


*2<Pi1)dx1dx? 


1  - 


rz 


-z  -  -z 


^Z^i' j«>dxidx2 


j 

rz  rz  r 

1  - 

J 

1 

N 

1 

N 

s  J 

rz 


;-[P..,P..,,P.., ]dx,x„dx_ 
y  ij  ij  JJ  12  3 


[Assuming  WLOG  that  j 


which  given  the  unknown  nature  of  is  the  simplest  form  this  can 

be  put  in. 

Theorem.  The  trivariate  method  gives  superior  results  than  does  the 
bivariate  method,  i.e.  (16)  when  applied  to  a  maximal  gives  a  lower 
upper  bou  i  for  PrfA]  than  does  (9)  applied  to  T  . 

Proof : 


(17)  l  Pr(B  )  =  I  Pr (A.  uA.)  =  l  Pr (A. )+Pr (A . )-Pr (A. nA. ) 

Ti  i.Jrto  1  J  1  J 

n 

=  y  Pr(A .)(#  edges  connected  to  A.)  -  V  Pr(A.nA.) 
i-1  1  1  i.jET0  1  J 


wlog  assume  j=i' 

(18)  7  Pr (B . . nB . , . , )  =  7  Pr(B..nB..,) 


ij  i  j 


ij  A  j  £T 


0 


ij.j'cT 


0 


L  Pr(A.)+Pr((A  uA. ,)  r*AT') 
i,  j  ,  j ,£Tq  J  ^11 


=  )  Pr(A.)(#  of  edges  connected  to  A.)~l 

i-1  1  1 


+  T  Pr(A.nA. ,)  n A? 


i  i  i ' cT  1  1  1 

1’1»1 


so  subtracting  (17)  from  (18)  shows  that  the  upper  bound  for  Pr(uAi) 


given  by  (16)  applied  to  is 


16 


[TPr (A .)*[//  edges  connected  to  A.]  -  Y  Pr(A.nA.)] 

1  i  4  r  T  ^  J 


ijcTr 


[EPr(A.)[/>  edges  connected  to  A.  on  T.-l]  +  Y  Pr  (A.  nA . , )  n  Ac ] 

1  0  jj'et;  1  J  J 


EPr (A  )  -  l  Pr(A.nA.)  -  5!  Pr((A.nA’.)  n  Ac) 

1  J  ijj’eT*  1  J  J 


EPr (A. )  -  Y  Pr (A. nA . ) 

*  «eTo  1  ' 


=  the  upper  bound  for  Pr(uA^)  given  by  (0)  applied  to  TQ. 


There  are  some  problems  with  using  (16)  on  to  find  an  upper 

bound  for  Pr[uA.]: 

l 

1)  It  may  not  be  practical  except  when  Ai,...,An  form  a  string  to 

find  a  best  tree  Tq.  Also  the  and  Tq  formed  above  do  not 

necessarily  give  the  lowest  upper  bounds  possible  for  Pr(uA_.)  when 
using  trees  with  bivariate  and  trivariate  normal  distributions . 

2)  It  is  computationally  very  expensive  and  time  consuming  to  calculate 

the  triple  interval  fZ  fZ  fZ  ^P,,,?  ,P.Jdx.x„dx,  using  the  IMSL 

“*Z  -z  — Z  1J  J  K  X  4  J 

procedure  DMLIN.  It  may  be  possible,  however,  that  superior  algorithms 
exist  for  integrating  trivariate  normal  distributions. 

3)  When  the  covariance  matrix  E  is  not  full  rank  there  is  no  computer 
algorithm  which  can  calculate  the  above  integral  (although  Uusipaika’s 
method  will  work  if  developed  into  a  computerized  procedure) . 

4)  Due  to  the  expensiveness  of  calculating  the  trivariate  normal 
intervals  it  would  be  impractical  to  use  the  Newton-Raphson  procedure 
with  this  method  to  find  Zq  for  a  given  otp-  Perhaps  some  other  type 
of  numerical  estimation  could  be  used,  see  chapter  V  part  b. 


Evaluating  the  Improvements  in  Upper  Bounds  for  and  U 


Estimates  for  z _  Given  by  Using  These  Methods 


It  is  envisioned  that  these  methods  could  be  applied  for  simul¬ 
taneous  inference  on  as  many  as  twenty  parameters  for  which  we  have 
jointly  normal  estimates.  It  has  been  shown  that  these  methods  give 
closer  upper  bounds  than  does  the  standard  Bonferroni  method.  The 
standard  Bonferroni  method  is  not,  however,  the  best  method  available 
for  calculating  upper  bounds  for  multivariate  normal  probabilities 
of  lying  in  cubes  centered  at  the  origin.  Sidak  proved  that 
if  x  is  an  n*n  vector  distributed  N((^),  Z) ,  where  Z  is  any 
n  *  n  covariance  matrix  with  diagonal  elements  of  1  then 


n  n 

Pr[  U  x  j£(-z,z)]  1  Pr[  U  y  i  (-z,z)] 
i-1  i-1 


where  y  is  an  (n * n)  vector  distributed  N((q),  I).  Therefore,  by 
assuming  the  components  of  any  vector  are  independent  and  calculating 
the  probability  that  the  vector  lies  outside  of  the  n  cube  centered 
at  the  origin,  one  gets  an  upper  bound  for  the  true  probability  that 
the  vector  lies  outside  the  cube.  This  upper  bound  is  smaller  than 
the  standard  Bonferroni  upper  bound  since  if  we  let 


A.  =  Pr{x.  i  (-z,z)};  A  =  u  {A.} 
11  i-1  1 


C.  =  Pr{y  i  (-z,z)};  C=  u{C} 
11  i=l 


Then  for  all  i:  Pr(A^)  =  Pr(C^)  and  Pr{C}  is  the  Sidak  upper  bound 
for  Pr(A)  while  ZPrCA^)  is  the  Bonferroni  upper  bound  for  Pr(A;  am 


< 


< 


Pr (A)  .  Pr{C}  ,  .IPr (C . ) 
Sidak  Bonferroni  1 


LPr (A^)  . 


This  conservative  independence  upper  bound  is  easily  calculated  since 
Pr{C}  =  l-(l-Pr{C^})n.  The  question  remains  as  to  when  the  methods  proposed 
in  this  report  give  lower  upper  bounds  than  does  the  conservative 
assumption  of  independence. 

The  answer  to  this  question  will  be  looked  at  in  two  different  ways. 

The  first  will  be  theoretically  to  see  if  we  can  determine  the  usefulness 

of  the  methods  as  relates  to  the  |P..ji,jeT_,  the  values  of  z  and  a 
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and  n  the  number  of  variables.  The  second  way  involves  empirical 
observations  of  improvement  for  particular  values  of  n,  ,  Zp  and  P 
where  P..  =  P^ ^  .ije^. 


a_.  Theoretical  Investigation  of  the  Amount  of  Improvement  in  Lowering 
a0  an(l  zq  Using  Improved  Bonferroni  Methods. 


Theorem  1.  For  any  tree  T^  and  any  collection  {P^..}  of  P„,ijeT 


ij 


ij 


0 


then  given  the  same  tree  Tq  and  the  same  {P„ }  except  replacing 

,  the  upper  bound  for 


5.  .  with  P!  .  where 

Vo  Vo 


Vo 


<  Ip, 


0J0 


Pr[uAp]  given  using  the  bivariate  and  possibly  the  trivariate  method  on 
the  new  set  of  {P..}  is  lower  than  the  upper  bound  for  Pr[uA.l 
given  using  the  respective  method  on  the  original  set  of  {P„}. 

Proof : 

(a)  Bivariate:  By  (Sidak)  (1  -  fZ  fZ  (P  )dx  dx~)  >  (1  - /Z  /z  4>,(P.'  .  )dx  dx  ) 

“  "Z  Z  Vo  -z  -z  z  iqj0  i  L 


iff  P. 


10J0 


<  ip! 


Vo 


Thus  the  upper  bound  given  by  (9)  will  be 


lower  for  the  new  set  of  {P_}. 


(b)  Trivariate  [may  not  always  be  true,  needs  to  be  nroven] .  Also 


need  to  consider  P...  where  A.--A.-A.-  on  the  tree. 

ij  i  i  j  • 

Theorem  2.  For  Jp  |  where  P^  <  1,  ijeT^.  Then  if  we  let 

A)  aQ(I)  be  the  a  calculated  for  a  fixed  Zp  under  assumption  of 

independence . 

ctg(B)  be  the  a  calculated  for  a  fixed  Zp  using  the  bivariate  method. 

a^CTr)  be  the  a  calculated  for  a  fixed  Zp  using  the  trivariate  method. 

Then  let  n  -*  00 .  There  exists  a  such  that  for  all  n  > 

ap (!)  1  !>f(B)  and  aF(I)  ^ap(Tr). 

Also,  if  we  let 

B)  2q(I)  be  the  z  calculated  for  a  fixed  af  under  assumption  of 

independence. 

Zq(B)  be  the  z  calculated  for  a  fixed  a^  using  the  bivariate 
method. 

Zg(Tr)  be  the  z  calculated  for  a  fixed  a^  using  the  trivariate 
method. 

Then  let  n  =  -  variables  -*•  *.  There  exists  a  Nq  such  for  all 
n  >  Nq:  Zp ( I )  <  Zp (B) ,  Zp(I)  <  Zp(Tr). 

Proof  of  A) . 

For  fixed  Zp  the  proof  is  easy  since  it  can  be  shown  that  for 
large  enough  n,  the  upper  bound  given  by  the  trivariate  method  is 
greater  than  1.  This  implies,  since  the  bivariate  method  is  never 
superior  to  the  trivariate  method,  that  the  bivariate  upper  bound  is 
also  greater  than  1. 


m 


i 


Formula  (19)  gives  the  trivariate  upper  bound  as 


IPr(A.) 


7  Pr  (A .  n  A . )  -  ^  Pr  (A .  n  A,  n  AC) 

i,jeT0  1  J  ij.jkonT^  1  J 


so  adding  a  new  element  A^  to  the  tree  adds  the  terms 


Pr (A  )  -  Pr (A  a A  )  -  Pr(A  ~ A,  n Ac) 
n  n  tn  n  k  m 


'  Pr'(Ann  <Ak  -  A„'C) 


and  if 


<  Pn  <  1  then  P  ,  ,  P,  ,  P  ranges  over  a 
-  0  nk  km  nm  & 


' nk 1 ’ 1  km 1  ’  1  nm 1 

compact  set  hence  Pr[A  n(A,  uA  )'"]  which  is  a  function  of  P  ,  ,  P,  , 

n  k  m  nk  kn 

P  and  is  bounded  bv  zero  must  achieve  its  minimum  value  at  some 
nm 

point  in  this  set.  But  for  every  point  in  this  set  Pr(A  A  (A,  jA  )c)  >0 

n  km 

since  x  =  x.  or  x  =  x  is  impossible.  Therefore  Pr[A  A 
n  k  n  m  1  n 

Q 

(A^  u  A^)  ]  >  £.  >  0.  Thus  adding  each  new  element  to  the  tree  adds 
at  least  £  to  the  trivariate  uppe  bound.  So  if  we  add  Q  >  t 
elements  to  the  tree  we  will  derive  the  trivariate  upper  bound  above 
1  and  1  >  (independence  upper  bound). 


Proof  of  B) 

This  is  very  complicated  and  will  not  be  done. 

QF.D. 


Two  Part  Theorem  Regarding  the  Relationship  Between  the  Improved  Bonferroni 
Lower  Bounds  for  Areas  on  Normal  (z)n  Cubes  Centerd  at  the  Origin  and  the 


Value  of  z. 


*««Vl * t HlmH '  1*1* .  u*j.« jUMi*  u<  i 


ij  i 


Part  I)  As  z  4  0  for  a  given  n  >  3  and  ncnsingular  correlation  sub- 

f  j'9'j  'l 

matrix  (I)’  of  size  q  >  3,  ana  x  *v  M  :  ’(E)  ,  there  exists  a 

q  U0J  nJ 

(depending  on  n  and  £)  such  that  z  <  z^  implies  that  the 
upper  bounds  given  for  Pr[X^ji(-z ,  z)  any  i]  are  greater  than  1  for 
both  the  bivariate  and  the  trivariate  methods. 

Proof: 

It  has  beer,  shown  already  that  the  bivariate  upper  bound  can  be 
no  lower  than  the  trivariate  upper  bound  so  it  suffices  to  prove  the 
statement  for  the  trivariate  upper  bound  onlv. 

Let  H  and  K  be  two  events  such  that  H  is  not  linked  to  K 
on  Tq ;  a  link  Bhj~Bjk  does  not  exist  on  for  any  j;  and  there 

exist  variables  Xn,XQ  such  that  [X  has  a  nonsingular 

distribution.  Then  the  trivariate  upper  bound  for  the  Pr[X  ^(-z,z)  any  i] 
will  equal 


Pr{  U  X.^(-z,z)i=l, . . . ,n} 
i=l  1 


+  Prt  n  (X  c(-z,z))n  X,  a  X ,£(-z,z) 

i/h,k  1  h  k 

+  Probabilities  of  other,  disjoint  terms 


>  [Pr{  n  (X.e(-z,z))nU(-zf2)nU(-z,z)!]  + 
i^h,k  1  n  * 


[1  -  Pr {  f)  X.  c(-z,z)  }] 
i-1 


=  l+[Pr{  f)  (X.c(~’,z))n  X,  t(-z,z)  r  x  tf(-z,z)} 
i*h,k  1  a  ^ 


Pr{  f\  X  e (-z , z) }  ] 
i=l  1 


But  as  z  *  0  the  ratio 


(Pri  P)  (X^e(-z,z))n  X^i^-z.z)  c  X^^(-z,z)}+  Pr{  A  X^e(-z,z)  }] 
i#h,k  i=l 


Pr  {  P\  X .  e  (-z  ,  z)  } 
i=l  X 


’  ^  X  e(-z,z); 
i^h.k 


Pr  ■.  ■'"N  X  .  c  (-z ,  z)  } 
i=l  1 


density  on  the  h,k  plane  passing  through  origin 
density  on  origin 


a 


=  *  due  to  nonsingularitv  of  distribution  of  (X  ,X„,X,  ,X,  ). 

'  p  U  n  k 


Therefore  there  exists  a  value  of  z  such  that  z  <  z  => 

Ju 


Pr {  H  X.£(-z,z)} 
i^h  tk  1 


Pr{  P.Xe(-z.z)} 
i-1 


1+Pr{  0  (X.e(-z,z))nx  *(-z,z)nX  ^(-z,z)  }-Pr{  f\x.£(-z,z)}  >  1 
i^h,k  K  -•-!  1 


trivariate  upper  bound  for  PrCz^C-z.z)  any  i}  >  1 


Part  II)  As  z  +  ®.  For  a  given  n  and  nondegenerate  (i.e.,  no  two 

events  are  the  same)  correlation  matrix  I  and  X  %  N(9),[I]  )  there 

0  n 

exists  an  z  depending  cn  n  and  £  such  that  z  >  z  implies  that 


Q  is  Che  upper  bound  estimate  of  the  probability  (X)  is  not  in 
the  cube  centered  at  zero  with  edge  length  z  using  independence. 
Q2  is  the  upper  bound  estimate  of  the  probability  (X)  is  not  in 
the  cube  centered  at  zero  with  edge  length  z  using  either  the 
bivariate  or  the  trivariate  method. 


<  (1  +  e)  . 


Proof .  [For  bivariate  method]. 

1)  First  we  will  show  that  for  a  given  tree  of  size  n 


Pr(A.r.A)  1 

lim  zt®  sup  • — - , — H—  <  -z—  for  all  z  >  0 
...  in 

i 


xrhere 


A.  =  Pr {X.jf  (-z,z)} 
A.  =  Pr-Xjrf  (-z,z)>. 


One  way  to  do  this  is  to  note  that  if 


f  \  r  f  ^  f  U 

Xi  0  |  Remember  [P..(^l  since 

'v  N  Q  |  J  |  ij 

X  U  I  P  1  £  is  nondegenerate. 

v  y  ^  '  ij  n 


WLOG  assume  P..  >0. 

ij  ~  _ 

Then  one  can  say  X.  =  P..X.  +  Fl-P^.(v)  where  v  is  N(0,1) 

J  ij  i  ij 

and  independent  of  X..  Now  for  all  c  >  0  and  o  >  0,  -•  z,  such 

l  -  b 

that  for  any  z  >  z. 


$ 


Pr{jX.j  >  z+a|jXil  >  z)  <  (~)  . 


So  for  anv  z  >  z. 


Pr{ |X.  >  z I  X.  >  z> 
3  1  1  1  i1 


£  Pr{|x.|  >  z+c] |X.|  >  z}  +  Pr{|X.|e[z>z+o]| |x.|  >  z} 
<_  P r {  ! X ^  |  >  z+o  |  jx^ |  >  z}  +  Pr{  X^  j z ,  z+o ]  ’  ,’x^|  >  z } 


£  I  <-1"Pii)  1  „  f  } 

-  ^+pr  y<— LzWpr^y>— J=—  f 
1  J  X-P7 .  (  il-P2.  J 

ij  13 
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aut  -  z  such  that  z  >  z  =*■ 
-  c  c 


f  -u+p„)  1 

Pr  <y  < - -i-  z> 


(*** ) 


1-P .  . 

Pr  Jv  >  - z  - 


1-P2. 

13 


ii  <k 


2  2  1 
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So  if  z  >  max(z,  ,z  )  then  (*),  (**)  and  (***) 
b  c 


■Vi  -v- .  v.>vv.A‘'.  <*.■ -*«yv  v.  ^  j 


pr{|x  j  >  z|  |X  |  >  2}  ^  ^ 

J  i  —on  —  2n 


=*  Pr  { lXj  |  >  zfl]x.|  >  z}  <  j-  pr{  |x  1  >  z} 

II)  So  WLOG  let  P„  =  sup^Jp^l  for  all  k.^eT^  and  z  be  larger  than 
max(zb,zc).  Then 

PrilxJ  >zH|X,|  >  z}  "  Pr{|x  i  >  z} 
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Q2  =  l  Pr{!XJ  1  z}  "  I  Pr{jx,|  >znix  I  >  z} 
k=l  k.leT .  K  1 
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1  Pr{  ;Xk'  >  z}  >_  Q 


L  Pr  i  U  ’Xk|  >  z} 

Independence  k=l 
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(1+e)  for 


The  proof  for  the  trivariate  method  follows  from  the  proof  of  the 
bivariate  method,  since: 


(I)  We  saw  in  the  bivariate  method 


;  Pr(A.nA.) 

lim  sup  -"p7jxt"-^  f°r  311  £<0, 

ZA°°i_l?j  1 
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which  implies  that 


Pr(A.nACnA'.)l 

,  .  1  i  i  t  e  .  , ,  „ 

lim  sup  - = — 7*\ —  ,  <  for  ail  c  >  0 

.  i  Pr  (A . )  -  2n 


(II)  So  the  trivariate  upper  bound  (Q^)  given  by  (19)  is 


(*)  XPr  (A  )  -  7  Pr(A.nA. )  -  7  Pr  ( (A.  r>A.  ,  )nAC)  . 

ijd0  1  3  ijj7sT£  1  3  3 


but  from  the  proof  of  the  bivariate  method  and  the  fact  that 
Pr((A^'Al )oa?)  _<  Pr(A^nA^)  it  follows  that  for  z  large  enough 


(*)  >  nPr(A.)-(n-l)  ~  Pr (A.)-(n-2)  ^  Pr(A.) 
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=  nPR(A  )  -  Pr(A. )  <  (n-c)Pr(A.) 
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=  U-“)  I  Pr( !  Xfc  i  >  z} 
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(III)  Now 


7  {PrjXjJ  >  z)  >  Q  >  Pr  (  u  | |  >  z) 

Tq  independence  k=l 

Qi  VPr,xJ >  z}  r  i  i 

q:  < - ol - <  ITT  " (1+£)  for  n  >  2  • 


Conclusion. 


The  previous  theorems  of  this  and  other  sections  give  us 
(a)  The  trivariate  method  always  gives  lower  upper  bounds  than  does 
the  bivariate  method. 


(b)  Everything  else  remaining  the  same  (i.e.  n,  T  ,1)  replacing  a 
Pi  i  on  the  tree  with  a  (P?2?  ]  such  that  Ipf1^  |  <  IpV  | 

Vo  Vo  Vo  1  Vo 

will  always  produce  a  lower  upper  bound  using  the  bivariate 
method  and  will  probably  also  produce  a  lower  upper  bound  using 
the  trivariate  method. 

(c)  For  fixed  n,  nonsignular  l ,  TQ  and  T'  making  zQ  close  to 
zero  causes  the  bivariate  and  trivariate  methods  to  produce  higher 
upper  bounds  in  relation  to  those  produced  by  the  assumption  of 
independence. 

(d)  For  fixed  n,  nonsingular  I,  TQ  and  T'  making  zQ  verv  large 
causes  the  bivariate  and  trivariate  methods  to  produce  upper  bounds 
which  get  infinitely  closer  to  those  produced  by  the  conservative 
assumption  of  independence. 

(e)  For  bounded  }P  |  ^ ;  increasing  n  eventually  causes  the 

bivariate  and  trivariate  methods  to  produce  poorer  upper  bounds 
with  respect  to  those  produced  by  the  conservative  assumption  of 
independence. 

Thus  we  can  see  there  is  a  complex  relation  involving  n,  (Zp  or  ip), 
and  the  P„  1  i,j  n  which  governs  the  effectiveness  of  using 
the  trivariate/bivariate  methods  over  the  assumption  of  independence. 

In  the  next  section  we  try  to  develop  an  actual  picture  of  this 
effectiveness  using  sample  calculations  of  upper  bounds  for  atn  given 


z  and  z  given  a  for  selected  values  of  n,  and  P..i,jeTn. 
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b.  Actual  Upper  Bounds  Obtained  for  ZqIo^  and  aQ 1 zp  From  Using 
the  Trivariate  and  Bivariate  Methods. 

1.  Table  1  contains  upper  and  exact  bounds  for  a,J  z^  obtained  bv 

u  F 

using  several  methods  involving  several  assumptions  on  Zp,  n,  and  I: 

(i)  Zp  was  taken  to  be  1.96,  2.25,  2.5  and 
felt  that  these  values  of  z  were  representative  of  what  would  be  of 
interest  in  practice. 

(ii)  n  was  taken  to  be  4, 5, 6,... 20  because  it  was  felt  that  in 
practice,  simultaneous  inference  on  4  to  20  variables  would  be  desired, 

(iii)  The  methods  used  were 

A.  Column  3:  Conservative  assumption  of  independence  [Corollary  3, 

Sidak,  1971]  which  says  for  X  •v  N(0,£) 
n  n 

Pr{  y  xV(-zFai,zFai)}  <  1  -  n  Pr{x±e (-ZyCK, ZpCn)} 
i=l  i=l 


B.  Column  4:  Bivariate  method  assuming  the  |P  |  all  have  the 

XiXj 

same  value  P  for  e^eT.  The  value  P  was  allowed 
to  be  .1, .2, . . . , .9, .99.  Calculations  were  done  by  a 
computer  program  using  formula  (9)  and  the  IMSL 
procedures  MDNOR  and  MDBNOR  to  calculate  bivariate  and 
univariate  normal  cube  probabilities  respectively 

C.  Column  5:  Trivariate  method  making  the  same  assumptions  as  in  B 

but  also  assuming  P  =  (P)  for  e...  ,eT’.  Calcula- 

x£xj»  ij’lj 

tions  were  done  by  computer  programs  using  formula  (16) 
and  the  IMSL  procedures  MDBNOR  and  DMLIN  to  calculate 
trivariate  and  bivariate  cube  probabilities  respectively. 
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Column  6:  The  exact  value  for  the  equicorrelation  case  where 

Pio  =  P  for  all  k,£..  The  formula  taken  from  Hahn  and 


Henderson  is: 


zF+*p  y 


-z„+/p  y 


Pr{  \J  A. 
i=l  1 


LJA.}  =1-  ff  /I^ 


-p  tf>(x)dx  - 


/l-p  P(x)dx  P(y)dy 


Again  P  was  allowed  to  be  . 1, .2, . . . , .9, .99 

E.  Column  7:  Trivariate  method  making  the  same  assumptions  as  in 

2 

B,  but  also  assuming  P  =  (P)  for  e..  ..,eT’  . 

x^Xj»  f  j  >  j  j 

This  is  an  AR(1)  time  series  process.  Calculations 
were  done  as  in  Column  5. 

Examination  of  this  table  will  give  some  ideas  as  to  the  improve¬ 
ment  of  the  upper  bound  for  <Xq j Zp  which  are  obtained  through  the 
bivariate  and  trivariate  methods.  The  two  points  of  interest  are: 

1.  How  much  lower,  if  any,  are  the  upper  bounds  for  obtained  from 
the  bivariate  and  trivariate  estimates  than  are  those  obtained  from  the 
conservative  assumption  of  independence. 

2.  How  close  are  the  bivariate  and  trivariate  upper  bound 
estimates  to  the  actual  values. 

Column  6  in  Table  1  gives  actual  values  of  for  the  case 

where  there  is  equicorrelation  among  all  the  variables,  i.e. 
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Column  3  gives  the  independence  upper  bound  for  the  in  column  6 

while  column  4  is  the  bivariate  upper  bound  for  the  in  column  6. 

Column  5  is  the  trivariate  upper  bound  for  the  in  column  6. 

The  pattern  observed  in  these  tables  is  far  too  complex  to  be 
described  in  the  limited  space  of  this  section.  But  some  general 
observations  should  be  made. 

a.  For  correlations  of  less  than  .3,  the  bivariate  and  trivariate 
methods  do  not  improve  much  or  at  all  over  the  conservative  assumption 
of  independence  for  all  values  of  n  and  Zp  examined. 

b.  For  P. .  >  .6,  z_  >  2.25  and  all  n,  the  bivariate  upper  bounds 

lj  —  F  — 

for  Oq  are  much  smaller  than  are  the  independence  upper  bounds.  Example 

z  =  2.75,  P..  =  .8,  n  =  18  the  independence  upper  bound  for  ct_  is 
F  lj  U 

.10201,  while  the  bivariate  upper  bound  for  is  .07575. 

c.  For  all  P  . ,  z  >  2.25  and  all  n  the  trivariate  upper  bound 

ij  t  — 

for  the  AR(1)  process  is  not  substantially  smaller  than  is  the  bivariate 
upper  bound.  Example  Zp  =  2.25,  P„  =  .4,  P  ,  =  .lt>  and  n  =  6.  The 
bivariate  upper  bound  for  is  .13502  and  the  trivariate  upper  bound 

for  is  .13254. 

d.  For  P  .  >  .4  and  all  z_  and  n,  the  trivariate  upper  bound 

i3  -  F 

for  in  the  equicorrelated  case  is  significantly  lower  than  is  the 

bivariate  upper  bound  for  yet  is  much  larger  than  the  actual  value 

of  a„  especially  if  P..  is  large.  Example  1:  P. .  =  .5,  n  =  8,  z  =  2.5. 


KM 


The  bivariate  upper  bound  is  .08998,  the  trivariate  upper  bound  is  .08400 


and  the  actual  value  is  .07456.  Example  2:  P..  =  .7,  n=13,  z_  =  2.25. 

ij  F 

The  bivariate  upper  bound  for  aQ  is  .23582,  the  trivariate  upper  bound 
is  .19754  and  the  actual  value  is  .13471. 

For  the  high  correlations,  i.e.  _>  .7  the  trivariate  and  the 

bivariate  methods  give  better  upper  bounds  for  for  n  even  as 

large  as  20  than  does  the  conservative  assumption  of  independence.  In 
the  case  where  the  correlations  between  variables  become  multiplicatively 
smaller  the  farther  away  one  moved  along  the  tree  (i.e.  AR(1))  the  reductions 
in  the  upper  bounds  from  using  the  trivariate  method  were  not  much  greater 
than  were  those  from  using  the  bivariate  method.  If  the  correlations 
of  variables  far  apart  from  each  other  on  the  tree  are  high  with  respect 
to  the  correlations  between  adjacent  variables  on  the  tree  (equi-correlated 
case)  then  even  the  trivariate  upper  bounds  for  were  much  higher 

than  the  true  values  of  cx^. 


2.  Tables  2A,  2B  and  2C  contain  upper  bounds  for  Zq  obtained  by 


using  the  bivariate  method  along  with  several  assumptions  on  n  and  2. 

(i)  3  was  taken  to  be  .05,  .10,  and  .20  because  it  was  felt 

r 

that  in  practice  jointly  simultaneous  confidence  intervals  of  80%, 

90%  and  95%  probability  would  be  desired. 

(ii)  n  was  taken  to  be  4, 5, 6,..., 20  because  it  was  felt  in  practice 
that  simultaneous  inference  on  4  to  20  variables  would  be  desired. 


(iii)  (2)  was  chosen  so  that 


p 

‘  X.X. 
1  1 


e„eT  always  equaled  some 


constant  P.  P  was  set  at  . 1 , . 2 , . . . , . 9 , .95  and  .99. 
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Zq | was  calculated  by  choosing  an  initial  Zg  and  then 

calculating  by  formula  (9)  using  the  IMSL  procedures  MDNOR  and 

Z0 

MBDNOR.  A  new  Zq  was  calculated  by  (9.5).  This  was  repeated  until 

la.  -a_ |  <  .0001  at  which  point  z„  was  set  to  z  . 

£  F '  0  0 

Z0 

No  attempt  was  made  in  this  table  to  find  trivariate  upper  bounds  for 
Zp!dg.  Due  to  the  present  computational  difficulty  of  working  with  tri¬ 
variate  intervals,  it  would  be  difficult  to  use  the  Newton-Raphson  procedure 
to  do  this  as  could  be  done  for  bivariate  upper  bounds.  Some  type  of 
interpolation  method  would  probably  be  the  best  way  to  obtain  these 
bounds.  One  crude  method  would  be  by  finding  Zg  (the  bivariate 
upper  bound  for  Zg,  Remember  zg  >  z^.  Calculate  (under  the  trivariate 


method)  |  Zg.  It  follows  that  Og  <  a^. 


Find  z  such  that  the 
u 


single  interval  probability  that  X^. l(-z^,z^)  is  an.  Then  calculate 


O' 


a  z 
u '  u 


for  all  of  the  intervals  together  using  the  trivariate  method.  It 


will  follow  that  a  >  a„.  Finally  let 

u  0 


(19) 


a  -a„  a  -n_ 

*  u  0  ,  0  B  „ 

Z„  =  -  Zu  +  -  Z 

F  a  -a„  B  a  -a„  u 


u  B 


u  B 


The  z^  we  obtain  will  be  smaller  than  that  obtained  bv  the  bivariate 
r 

method.  It  also  should  be  larger  than  the  true  trivariate  upper  bound 

z„  since  it  makes  sense  to  believe  that  z 
F  a 


is  a  convex  function,  i.e. 


As  was  the  case  in  Table  1;  the  patterns  shown  in  Tables  2A,  2B 
and  2C  are  too  complex  to  describe  completely.  Some  important  things 
can  and  should  be  noted  from  these  tables,  however.  First  for 
0^  =  .05,  .10  and  ,20;  the  bivariate  method  generally  gives  lower 
values  of  z^  for  all  n  than  does  the  conservative  assumption  of 
independence.  Although  the  reduction  in  values  may  seem  small  they 
are  actually  larger  when  one  considers  the  reduction  in  the  hyper¬ 
volume  of  the  cube  as  opposed  to  the  reduction  in  the  edge  length  of 
the  cube. 

Example  1:  a_  =  .20,  n  =  6 ,  P . .  =  .4. 

F  lj 

The  conservative  assumption  of  independence  gives  an  upper  bound  for 
Zq  of  2.111  while  the  bivariate  method  gives  an  upper  bound  for  Zg 
of  2.087.  This  representates  a  reduction  of  edge  length  to  .9886  of 
its  original  length.  But  it  also  represents  a  reduction  in  the  hyper 

«Hf 

volume  of  the  cube  to  (.9886)  =  .9337  of  its  original  value. 
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Example  2:  =  .20,  n=10,  P  =  .7  . 

The  conservative  assumption  of  independence  gives  an  upper  bound  for 
Zq  of  2.309  while  the  bivariate  method  gives  an  upper  bound  for  z^ 
of  2.213.  This  represents  a  reduction  in  edge  length  to  (.9584)  of 
its  original  length.  But  it  represents  a  reduction  in  the  hyper  volume 
of  the  cube  to  (.9584)  ^  =  .6540  of  its  original  volume. 

Example  3:  =  .05,  n=15,  =  .9  . 

The  conservative  assumption  of  independence  gives  an  upper  bound  for 
Zq  of  2.932  while  the  bivariate  method  gives  an  upper  bound  for  z^ 
of  2.742.  This  represents  a  reduction  of  edge  length  to  (.9352)  of 
its  original  length  and  a  reduction  in  hyper  volume  of  the  cube  to 
(.9352)  3  =  (.3660)  of  its  original  volume 

Example  4:  a^.  =  .10,  n=20,  P  =  .99  . 

The  conservative  assumption  of  independence  gives  an  upper  bound  for 

Zq  or  2.799  while  the  bivariate  method  gives  an  upper  bound  for  z^ 

of  2.216  which  is  .7917  of  its  original  length  and  a  reduction  in  hyper 

**20 

cube  column  to  (.7917)  =  .00936  of  its  original  volume. 

Table  3  compares  trivariate  and  bivariate  upper  bounds  for 

Zfja  =  .10  for  n  =  5,10  and  20  and  P..  =  .6, .7, .8  and  .9.  The 
Ur  i  j 

bivariate  upper  bounds  are  from  Table  2B.  The  tri.variate  upper  bounds 

are  calculated  as  in  (19)  using  both  the  case  w.  ere  P..'  =  (P..) 

13  il 

2 

e..  ..,eT'  and  the  case  where  P  ,=  (P..)  e..  ..,eT'.  The  tri¬ 
ll, 33  ij'  13  13  >33 

variate  method  does  give  significant  reductions  below  the  bivariate 

method  in  the  case  where  P . . ,  =  P . . . 
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Example  1:  a  =  .10,  n=  5,  P..  =  .7,  P..,=  .7. 

F  13  JJ 


The  trivariate  upper  bound  for  z^  is  2.17  while  the  bivariate  upper 


bound  for  z^  is  2.23.  This  gives  an  edge  length  ratio  reduction  of 


.973  and  a  hypercube  volume  reduction  of  (.973) 


**5 


.872. 


Example  2:  a  =  .10,  n  =  20,  P..  =  .9,  P  =  .9. 

F  13  j3 


The  bivariate  upper  bound  for  z^  is  2.59  while  the  trivariate  upper 
bound  for  Zq  is  2.49.  This  gives  an  edge  length  ratio  of 


**20 


.961  and  a  hyper  volume  reduction  of  (.961)  =  .455. 


In  the  case  where  Fjj'  =  (Fij)  the  reductions  given  by  the 


trivariate  method  are  not  as  impressive. 


Example  3:  a  =  .10,  n  =  10,  P . .  =  . 8,  P . , ,  =  .64. 
F  13  33 


The  bivariate  upper  bound  for  z^  is  2.44  while  the  trivariate  upper 
bound  for  z^  is  2.42.  This  gives  an  edge  length  ratio  of 
.992  and  a  hyper  volume  reduction  of  (.992)  ^  =  .921. 


C.  Summary 


In  part  a.  of  this  section  it  was  shown  that  the  bivariate/ 
trivariate  methods  gave  better  upper  bounds  with  respect  to  those 
given  by  the  conservative  assumption  of  independence  as  absolute  values  of 


correlations  increased,  number  of  decreased  and  z  became  neither 


too  large  nor  too  close  to  zero.  In  part  b  it  was  shown  that  for  large 


absolute  correlations  between  X..,  e..cT  the  bivariate  method  gave 

13  13  6 


considerably  lower  upper  bounds  for  a  I z  and  zA j a  than  did 

Ur  Ur 


independence.  It  was  also  shown  that  if  P..i  was  large  for  e..  .  .,cT' 

1J  13.3J 


that  the  trivariate  method  gave  considerably  lower  upper  bounds  for 


36 


was 


a0  z  and  zn  a_  than  did  the  bivariate  method.  If  P  , 
u 1  F  0 1  F  i  j ' 


small  compared  to  P..,  i.e.  P..,=(P..)  then  the  upper 

ij  '  ij’  ij 


bounds  for  ao'ZF  an<*  ZF^a0  8^ven  by  the  trivariate  method,  were 
not  much  smaller  than  those  given  by  the  bivariate  method.  It 
is  felt  that  the  improvements  in  the  upper  bounds  were  enough  to 
warrant  usage  of  these  methods  in  practice. 


It  is  desired  to  obtain  confidence  intervals  for  B.  where 

1 


i  =  1 , . . . ,n.  The  current  method  used  to  find  the  largest  z^  such 
that 


(21) 


Pr 


u  B.  i  (B.-zn(X’X)‘^o,  B.+zn(X'X)Tjo) 

J^=-^  1  1  U  11  1  U  11 


^aF 


is  the  standard  Bonferroni  method  which  involves  finding  z^  such  that 


Pr[B.  i  (B.-z_.  (X’X) .  fa,  B.+zn(X'X) . fa) ]  =  a  /n.  Smaller  values  of 
1  1  U  11  1  U  il  r 


Zq  which  will  have  property  (21)  can  be  obtained  by  using  (8)  to 


form  the  "Best  Bivariate  Tree:  Tq"  then  iteratively  using  the  Newton- 


Raphson  procedure  (9.5)  with  f(z)  =  (9)-ap  and  f'(z^)  =  (10) 


Even  smaller  values  for  z^  could  be  obtained  by  applying  (16)  in 


some  fashion  to  the  tree  Tq  created  from  Tq  by  (13) . 


It  is  important  to  note  that  often  o  is  not  known  and  is 

~  2 
(Y . -Y) 


estimated  by 


(n-1) 


which  then  means  that  confidence  intervals 


for  B.  are  created  using  a  instead  of  a.  All  of  the  methods 
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described  previously  could  then  be  easily  modified  to  use  multivariate 
t  distributions  instead  of  multivariate  normal  distributions. 

b .  Time  Series. 

These  simultaneous  inference  procedures  can  be  applied  to  estima¬ 
tions  of  probabilities  of  events  involving  autoregressive,  moving 
average  and  other  types  of  time  series  models.  For  instance  let's 
look  at  the  MA(3)  case  where 

Yi  =  gi(Xi)+82(Xi-l)+g3(Xi-2)  where  1  =  time;  8j.g9.g3 

2  2 

are  known  constants,  and  the  are  iid  N(0,o  ) ,a 

known . 

The  goal  is  to  produce  simultaneous  at  least  (l-aF)  confidence 
intervals  for  the  values  of  [Y^,Y^+^, .  . . ,  Yi+n_j  1  given  no  prior 
knowledge  of  X^,...,X^  or  Y^,...,Y^  Without  loss  of  generality 
let  a2  =  1  and  [gj+g^+g^]  =  !♦  Then 

where 

p2  =  f  gig3 ] 


n 

Using  the  notation  A  =  [  u  X.  t  (~z,z)]  and  [A.  =X.  i  (-z,z)]  then 

i=l  1  11 

the  best  bivariate  tree  Tn  is  [A. -A. , • *-A. ,  ].  The  standard 

0  1  l+l  l+q 

ct 

Bonferroni  z^  is  the  largest  z  such  that  Pr (X^£ (-Zq,Zq)  £  “• 

Smaller  values  for  z^  could  be  computed  iteratively  by  using  the  Newton- 
Raphson  procedure  (9.5)  on  T^  with  f  =  (9)-oip  and  f'  =  (10).  Still 
smaller  values  for  z^  could  be  computer  by  applying  (16)  in  some  fashion 

to  the  tree  T^  created  from  Tq  by  (13). 


Even  smaller  values  for  z could  be  obtained  by  applying  (16)  in  some 
fashion  to  the  tree  created  from  by  (13).  In  this  case  we 

might  need  Uusipaika's  formula  to  integrate  over  singular  three-dimensional 
distributions  since  we  are  dealing  with  linear  combinations  which  could 
mean  singularity  in  the  covariance  matrices  of  three  adjacent  points 


along  Tq. 


d.  Survival  Times.  (A  special  application). 


Up  until  now  the  exact  correlation  structure  of  x^,x^,  i.jcTg 

has  been  assumed  to  be  known.  Below  is  presented  an  application  of 

the  bivariate  method  in  a  situation  where  the  exact  correlations  are 

not  known.  The  application  will  be  quite  complex  and  given  without 

rigorous  justification.  The  author  plans  to  give  more  justification 

in  a  later  technical  report  or  thesis.  Those  who  are  not  interested 

in  this  application  may  wish  to  skip  this  section. 

Assume  the  classical  survival  model  with  no  censoring  where  the 

time  of  death  is  a  random  variable.  Let  S  be  the  probability  of 

i 

surviving  to  time  t^  where  0  <  t^  <  t^  <  • • •  <  t^.  Assume  that  a 

sample  of  m  individuals  is  followed  starting  at  time  zero  and  the 

time  of  death  of  each  individual  is  recorded.  Let  S  =  (number  of 

1 

individuals  alive  at  time  t .) /m  be  the  sample  estimate  of  S 

1  i 

for  i  =  l,2,...,k.  If  m  is  large  then  by  the  multivariate  central 
limit  theorem 


S  -S’ 
C1  C1 

s  -  S 

C2  tZ 


S  -  S 
.  Ck  CkJ 


approximately  N 


where  Z  =  S  (1-S  ) 

11  Ci  i 
and 


(1-St  )St  for  i  <  j  . 

i  j 
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Often  it  is  desirable  to  build  (l-a^)  joint  confidence  intervals 

for  the  ,  i=l,2,...,k.  Unfortunately,  the  variance  of  (S  -S  ) 
i  Ci  Ci 

and  the  correlation  of  (St  -St  ,  S£  -St  )  depend  on  the  true  unknown 

i  i  j  j 

values  of  S  and  S 
1 

The  current  way  to  construct  these  simultaneous  confidence 

intervals  uses  the  Bonferroni  method.  For  each  time  t^ ,  i=l,...,k, 

upper  and  lower  limits  for  the  confidence  interval  of  S  are 

i 

derived  in  (Fleiss  page  14)  and  given  below 


f 

ft: 

•ft 

$ 

^1 

S$ 

W; 

I 


(22a) 


(22b) 


where 


;  .  !i!!i 

’t  .5. 
u,  1  1 


B . — <t> 
_ i _ i 

Ci.,i  si 


(St  is  the  upper  limit) 


:  (St  is  the  lower  limit) 

'  ,  i 


B  =  2mS  +  (z  )2-l 

i  UF/2k' 


i  .  =  Z ,  . 

1  aF/2k 


<Z<“t/oi,)>2"<2  +  “)  +  «  <mU  -  S  )+l) 

F/2k  l  l 


H  =  2(m+(z  )2) 

k  F/2k; 


(aF/2k)  aF 

and  z,  ,  is  the  value  such  that  i}>(x)dx  =  (1-—) 

<-aF/2k; 

The  z,  N  term  in  the  above  formulas  was  derived  from  the 
Bonferroni  procedure:  dividing  k.  Below  is  a  procedure  based 

on  the  bivariate  method  which  derives  a  smaller  value  for  the  z  term 
producing  tighter  intervals  which  will  have  at  least  a  (l-a^)  joint 
confidence. 


1 

I 


STEP  I. 


t  <<fe  ' 


|g$j 


Clearly  the  best  tree  Tn  for  this  problem  is  S  -S  -...-S 

C1  .2  Ck 

The  sample  estimate  of  the  correlation  between  S  and  S  is 

i  i+1 


Cor (S  ,S  )  =  - 
i  i+1 


St  “^t.’ 
l+l  1 


/s  ( i— s  )  s  a-sr  ) 

7  i+1  i+l  l  i 


for  i=l,2 , . . . ,k-l 


So  using  these  correlation  estimates,  derive  zo!ap  from  the 

bivariate  method  applied  to  T^  using  the  Newton-Raphson  procedure 

with  f(z)  =  (9)-a  and  f'(z)  =  (10).  Replace  z,  .  by  z. 

F  (aF/2k)  0 

in  equations  (22a)  and  (22b)  and  let  the  k  upper  limits  obtained 

be  S'  ,S'  ,...,S'  and  the  k  lower  limits  obtained  by 

tlu  C2u  ku 

c  »  Cl  Cl 

J  , . . . 

li  21  ki 


STEP  II. 

(a)  Based  on  the  values  S'  ,S'  ,  ...,S'  calculate 

lu  C2u  ku 

cor(S  ,  S  )  using  formula  (23). Then  using  these  calculations, 

iu  ti+l,u 

derive  zolap  from  the  bivariate  method  applied  to  Tq  using  the 

Newton-Raphson  procedure  with  f(z)  =  (9)-aj,  and  f'(z)  =  (10). 

Replace  z,  .  by  this  value  of  zn  in  equation  (22a)  to  obtain 
'‘aF/2k;  U 

new  upper  limits  S"  ,S"  ,...,S" 

iu  C2u  Cku 


&*• 


P 

IW 

row 

[Ay 

§ 

m 


Kt? 

m 


y  •  •  •  > 


calculate 


(b)  Based  on  the  values  S  ,S 

tU  C2i 


~ki 


cor(s'  ,s'  )  using  formula  (23).  Then  using  these  correlations, 

U  *1+1,8. 

derive  zo!ap  from  the  bivariate  method  applied  to  using  the 

Newton-Raphson  procedure  with  f(z)  =  (9)-Op  and  f'(z)  =  (10). 

Replace  z,  >  by  this  value  of  z  in  equation  (22b)  to  obtain 

(aF/2kJ  U 

new  lower  limits  S"  ,  S"  ,  ...,S" 

Cl2,  t2l  tk5- 

It  is  claimed  that  the  (S'J  ,S"  )  will  form  a  set  of  joint 

tot. 

1U 

confidence  intervals  of  level  at  least  (1-cip)  and  that  these  intervals 

will  be  tighter  than  those  produced  using  the  Bonferroni  method. 

The  justification  for  this  claim  is  that  in  Step  1,  we  calculate 

the  upper  limits  S'  ,...,S'  and  use  these  upper  limits  to  get  a 

lu  Cku 

better  estimate  of  the  true  upper  bound  correlation  structure  and  with 

this,  get  the  estimates  S”  ,...,S"  in  Step  II.  Likewise  for  lower 

lu  Cku 

limits. 
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VII.  Modifications  on  Regions,  Distributions  and  Methods . 

The  improved  Bonferroni  method  described  which  used  integration 

over  univariate,  bivariate  and  trivariate  distributions  to  get  upper 
bounds  for  probabilities  over  multivariate  normal  cubes  centered  at 
zero  could  be  slightly  modified  to  handle  the  following  situations. 

A.  Everything  described  previously  could  be  applied  to  handle 


r#v;* 

I'tyi 

Mjf. 

•5v. 

& 


Mi 

W*'. 

frm 

mi 

i 


multivariate  t  distributions  by  replacing  the  univariate,  bivariate 
and  trivariate  normal  integrals  by  the  appropriate  univariate,  bivari¬ 
ate  and  trivariate  t  integrals.  The  conservative  assumption  of 
independence,  as  an  upper  bound  for  the  area  outside  the  rectangle 


II 

•V*.* 

VM 


will,  however,  not  necessarily  hold  for  t-distribution  regions.  Also 

the  correlation  coefficients  P,  „  must  be  known  (as  will  be  the 

k,  x, 

case  for  linear  regression  estimates  of  coefficients).  Further  studv 
needs  to  be  done  for  the  case  where  P^  ^  are  estimated. 

B.  The  cubes  centered  at  zero  could  be  replaced  with  rectangular 
regions  centered  anywhere.  Analogs  of  (9)  and  (16)  could  then  be 
used  replacing  the  lower  limit  of  integration  -z  with  z^  and  the 
upper  limit  of  integration  z  with  z^_,  for  the  i-th  marginal. 

Unless  the  rectangle  is  centered  exactly  at  the  origin,  the  conservative 
assumption  of  independence  as  an  upper  bound  for  the  area  outside  of  the 
rectangle  will  not  necessarily  hold. 

C.  Some  or  all  of  the  rectangular  regions  may  be  replaced  with 
half  plane  or  orthant  like  regions  in  which  case  analogs  of  (9)  and 
(16)  are  used  with  either  the  lower  limit  of  integration  -z  being 
replaced  by  -®  or  the  upper  limit  of  integration  z  being  replaced 
by  tx>.  Again  the  conservative  assumption  of  independence  as  an  upper 
bound  for  the  area  outside  of  the  region  will  not  necessarily  hold. 


i 

I 


s' 

SXv 


a 


D.  These  methods  could  be  modified  to  work  for  unions  of 
disjoint  sets  of  the  types  in  B  and  C  with  normal  or  t 
distributions.  The  calculations  might  be  quite  complicated  however. 
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DIRECTORY  OF  TABLES  1 


n 

Upper  and  Exact  Bounds  for  =  Pr  (J  X .  I  (-o^Zp  .o^p)  where  X^N(0,I) 

i=l 

for  zp  =  1.96,  2.25,  2.5,  2.75 
n  =  4,5,6, . . . ,20 
Calculated  under: 

Column  (3) :  Conservative  assumption  of  independence 

Column  (4):  Using  bivariate  method  assuming  there  exists  a  tree  T 

connecting  (X  ,X~,...,X  )  with  |P  !  =  P  for  all  e  ,-:T 
12  n  1  xixj 

Column  (5):  Using  trivariate  method  assuming  there  exists  a  tree  T 

as  described  for  column  4  and  a  tree  T'  constructed  from 

that  tree  such  that  P  =  P  for  all  e.^  . .,£T'. 

ij 

Column  (6):  The  exact  probability  calculated  under  the  assumption  of 
eauicorrelation  P  s  P  for  all  k;£ 

Vi 

Colmun  (7):  Using  trivariate  method  assuming  there  exists  a  tree  T 

as  described  in  column  4  and  a  tree  T'  constructed  from 

2 

that  tree  such  that  P  =  P  for  all  e..,e..,sT. 


i'LViViW 


Table  la.  Upper  Bounds  For  Pr((  (J  X.)<(-z,z)} 
'  ‘  i=l 


z=l . 96 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariatc 

P  =  p 

ij'  ij 

ij.jj'eT 

Equicorrelated 

Case  P.  =  P 
k£ 

all  k, £ 

P .  .  ,  =  (P  .  . 
ij  i: 

ij.jj'eT' 

.ino 

3 

0.14261 

_ 

0.14188 

0.14188 

0.14212 

.100 

4 

0.18548 

— 

— 

0.18412 

_ 

.100 

5 

0.22620 

— 

— 

0.22408 

— 

.100 

6 

0.2643S 

— 

— 

0.26190 

- — 

.100 

7 

0.30164 

— 

— 

0.29771 

_ 

.100 

8 

0.33656 

— 

— 

0.33165 

— 

.100 

9 

0.36973 

— 

— 

0.36382 

— 

.100 

10 

0.40124 

— 

— 

0.39432 

- — 

.100 

11 

0.43117 

— 

— 

0.42325 

— 

.100 

12 

0.45961 

— 

— 

0.45069 

— 

.100 

13 

0.48663 

— 

— 

0.47674 

— 

.100 

14 

0.51299 

— 

— 

0.50147 

— 

.100 

15 

0.53668 

— 

— 

0.52496 

— 

.100 

16 

0.55984 

— 

— 

0.54727 

— 

.100 

17 

0.58185 

— 

— 

0.56846 

- - 

.100 

18 

0.60275 

— 

— 

0.58861 

— 

.100 

19 

0.62261 

— 

— 

0.60775 

— 

.100 

20 

0.64148 

— 

— 

0.62597 

— 

.200 

3 

0.14261 

_ 

0.13972 

0.13972 

0.14059 

.200 

4 

0.18548 

— 

0.18300 

0.18047 

0.18475 

.200 

5 

0.22620 

— 

— 

0.21827 

— 

.200 

6 

0.26489 

— 

— 

0.25379 

— 

.200 

7 

0.30164 

— 

— 

0.28721 

— 

.200 

8 

0.33656 

— 

— 

0.31871 

— 

.200 

9 

0.36973 

— 

— 

0.34845 

— 

.200 

10 

0.40124 

— 

— 

0.37656 

— 

.200 

11 

0.43117 

— 

— 

0.40315 

— 

.200 

12 

0.45961 

— 

— 

0.42835 

— 

.200 

13 

0.48663 

— 

— 

0.45224 

— 

.200 

14 

0.51229 

— 

— 

0.47493 

— 

.200 

15 

0.53668 

— 

— 

0.49648 

— 

.200 

16 

0.55984 

— 

— 

0.51697 

— 

.200 

17 

0.58185 

— 

— 

0.53647 

— 

.200 

18 

0.60275 

— 

— 

0.55039 

— 

.200 

19 

0.62261 

— 

— 

0.57274 

— 

.200 

20 

0.64148 

— 

— 

0.58961 

— 

II 

means 

value  is  high 

er  than  that 

given  by  con 

servative  assumption  of 

independence . 


Table  la  Continued. 


z-1.96 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariat 

P,  -  P.  . 

Equicorrelated 

P.  .  ,  =  (P 

ij  ij 

ij.jj’eT 

Case  P,  „  =  P 

all  k,e 

ij .jj'eT 

.300 

3 

0.14261 

0.14023 

0.13623 

0.13623 

0.13802 

.300 

4 

0.18548 

0.18535 

0.17734 

0.17409 

0.18092 

.300 

5 

0.22620 

— 

0.21846 

0.20891 

0.22382 

.300 

6 

0.26489 

— 

0.25959 

0.24127 

— 

.300 

7 

0.30164 

— 

0.30069 

0.27146 

— 

.300 

8 

0.33656 

— 

— 

0.29971 

— 

.300 

9 

0.36973 

— 

— 

0.32624 

— 

.300 

10 

0.40124 

— 

— 

0.35122 

— 

.300 

11 

0.43117 

— 

— 

0.37479 

— 

.300 

12 

0.45961 

— 

— 

0.39708 

— 

.300 

13 

0.48663 

— 

— 

0.41821 

— 

.300 

14 

0.51229 

— 

— 

0.43827 

— 

.300 

15 

0.53668 

— 

— 

0.45733 

— 

.300 

16 

0.55984 

— 

— 

0.47549 

— 

.300 

17 

0.58185 

— 

— 

0.49280 

— 

.300 

18 

0.60275 

— 

— 

0.50932 

.300 

19 

0.62261 

— 

— 

0.52511 

— 

.300 

20 

0.64148 

— - 

~ — 

0.54021 

- - 

.400 

3 

0.14261 

0.13646 

0.13146 

0.13146 

0.13427 

.400 

4 

0.15848 

0.17970 

0.16968 

0.16582 

0.17531 

.400 

5 

0.22620 

0.22293 

0.20791 

0.19658 

0.21635 

.400 

6 

0.26489 

— 

0.24614 

0.22557 

0.25740 

.400 

7 

0.30164 

— 

0.28436 

0.25194 

0.29844 

.400 

8 

0.33656 

— 

0.32259 

0.27645 

— 

.400 

9 

0.36973 

— 

0.36082 

0.29932 

— 

.400 

10 

0.40124 

— 

0.39905 

0.32075 

— 

.400 

0.43117 

— 

— 

0.34092 

— 

.400 

12 

0.45961 

— 

— 

0.35995 

— 

.400 

13 

0.48663 

— 

— 

0.37796 

— 

.400 

14 

0.51229 

— 

— 

0.39505 

— 

.400 

15 

0.53668 

— 

— 

0.41130 

— 

.400 

16 

0.55984 

— 

— 

0.42678 

— 

.400 

17 

0.58185 

— 

— 

0.44156 

— 

.400 

18 

0.60275 

— 

— 

0.45568 

— 

.400 

19 

0.62261 

— 

— 

0.46921 

— 

.400 

20 

0.64148 

— 

— 

0.48218 

— 

ij' 


" - "  means  value  is  higher  than  that  given  by  conservative  assumption  of 

independence.  ^ 


Table  la  Continued 


z=1.96 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariate 

P.  .,  -  P.  . 
ij  1 3 

ij . J J  * 

Equicorrelated 

Case  P,  .  s  P 
k£ 

all  k,£ 

V  =  (pi 

ij.jj'eT’ 

.500 

3 

0.14261 

0.13348 

0.12543 

0.12543 

0.12922 

.500 

4 

0.18548 

0.17223 

0.16013 

0.15582 

0.16769 

.500 

5 

0.22620 

0.21297 

0.19482 

0.18285 

0.20617 

.500 

6 

0.26489 

0.25371 

0.22952 

0.20728 

0.24465 

.500 

7 

0.30164 

0.29445 

0.26421 

0.22959 

0.28313 

.500 

8 

0.33656 

0.33520 

0.29890 

0.25013 

0.32160 

.500 

9 

0.36973 

— 

0.33360 

0.26917 

0.36008 

.500 

10 

0.40124 

— 

0.36829 

0.28693 

0.39856 

.500 

11 

0.43117 

- — 

0.40299 

0.30356 

— 

.500 

12 

0.45961 

— 

0.43768 

0.31920 

— 

.500 

13 

0.48663 

— 

0.47238 

0.33397 

— 

.500 

14 

0.51229 

— 

0.50705 

0.34796 

— 

.500 

15 

0.53668 

— 

— 

0.36125 

— 

.500 

16 

0.55984 

— 

— 

0.37390 

— 

.500 

17 

0.58185 

— 

— - - 

0.38597 

— 

.500 

18 

0.60275 

— 

— 

0.39751 

— 

.500 

19 

0.62261 

— 

— 

0.40857 

- — 

.500 

20 

0.64148 

0.41918 

.600 

3 

0.14261 

0.12509 

0.11809 

0.11809 

0.12263 

.600 

4 

0.18548 

0.16263 

0.14864 

0.14404 

0.15771 

.600 

5 

0.22620 

0.20018 

0.17919 

0.16677 

0.19280 

.600 

6 

0.26489 

0.23773 

0.20974 

0.18681 

0.22788 

.600 

7 

0.30164 

0.27527 

0.24029 

0.20498 

C. 26297 

.600 

8 

0.33656 

0.31282 

0.27084 

0.22154 

0.29806 

.600 

9 

0.36973 

0.35037 

0.30139 

0.23678 

0.33314 

.600 

10 

0.40124 

0.38791 

0.33194 

0.25090 

0.36823 

.600 

11 

0.43117 

0.42546 

0.36249 

0.26406 

0.40331 

.600 

12 

0.45961 

— 

0.39303 

0.27639 

0.43840 

.600 

13 

0.48663 

— 

0.42358 

0.28798 

0.47349 

.600 

14 

0.51229 

— 

0.45413 

0.29893 

0.50957 

.600 

15 

0.53668 

— - 

0.48468 

0.30931 

— 

.600 

16 

0.55984 

— 

0.51523 

0.31917 

— 

.600 

17 

0.58185 

— 

0.54578 

0.32857 

— 

.600 

18 

0.60275 

— 

0.57633 

0.33754 

— 

.600 

19 

0.62261 

— 

0.60688 

0.34613 

— 

.600 

20 

0.64148 

— 

0.63743 

0.35437 

— 

— "  means  value  is  higher  than  that  given  by  conservative  assumption 
of  independence . 
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Table  la  Continued 


2=1.96 


RHO  Number  Independence  Bivariate 


Trivariate 
P  =  P 

ij'  U 
ij  » ji ' £T ' 


Exact  for 

Equicorrelated 

Case  P,  s  P 
k£ 

all  k,  £ 


Trivariate 

V  •  ‘V 

ij.Jj'cT' 


.700 

3 

0.14261 

0.11689 

0.10920 

0.10920 

0.11414 

.700 

4 

0.18458 

0.15034 

0.13496 

0.13035 

0.14483 

.700 

5 

0.22620 

0.18379 

0.16072 

0.14839 

0.17553 

.700 

6 

0.26489 

0.21724 

0.18648 

0.16416 

0.20623 

.700 

7 

0.30164 

0.25069 

0.21224 

0.17820 

0.23692 

.700 

8 

0.33656 

0.28414 

0.23799 

0.19087 

0.26762 

.700 

9 

0.36973 

0.31759 

0.26375 

0.20243 

0.29841 

.700 

10 

0.40124 

0.35104 

0.28951 

0.21306 

0.32901 

.700 

11 

0.43117 

0.38449 

0.31527 

0.22292 

0.35970 

.700 

12 

0.45961 

0.41794 

0.34103 

0.23211 

0.39040 

.700 

13 

0.48663 

0.45139 

0.36679 

0.24071 

0.42109 

.700 

14 

0.51229 

0.48484 

0.39254 

0.24881 

0.45179 

.700 

15 

0.53668 

0.51828 

0.41830 

0.25646 

0.48248 

.700 

16 

0.55984 

0.55173 

0.44406 

0.26371 

0.51318 

.700 

17 

0.58185 

— 

0.46982 

0.27059 

0.54387 

.700 

18 

0.60275 

— 

0.49558 

0.27716 

0.57457 

.700 

19 

0.62261 

— 

0.52134 

0.28343 

0.60526 

.700 

20 

0.64148 

0.54709 

0.28943 

0.63596 

.800 

3 

0.14261 

0.10611 

0.09822 

0.09822 

0.10307 

.800 

4 

0.18548 

0.13416 

0.11838 

0.11413 

0.12809 

.800 

5 

0.22620 

0.16222 

0.13855 

0.12734 

0.15310 

.800 

6 

0.26489 

0.19028 

0.15871 

0.13867 

0.17812 

.800 

7 

0.30164 

0.21833 

0.17888 

0.14861 

0.20314 

.800 

8 

0.33656 

0.24639 

0.19904 

0.15747 

0.22815 

.800 

9 

0.36973 

0.27444 

0.21921 

0.16548 

0.25317 

.800 

10 

0.40124 

0.30250 

0.23927 

0.17279 

0.27819 

.800 

11 

0.43117 

0.33056 

0.25954 

0.17952 

0.30320 

.800 

12 

0.45961 

0.35861 

0.27970 

0.18576 

0.32822 

.800 

13 

0.48663 

0.38667 

0.29987 

0.19167 

0.35324 

.800 

14 

0.51229 

0.41472 

0.32003 

0.19701 

0.37825 

.800 

15 

0.53668 

0.44278 

0.34020 

0.20213 

0.40327 

.800 

16 

0.55984 

0.47084 

0.36036 

0.20697 

0.42829 

.800 

17 

0.58185 

0.49889 

0.38053 

0.21155 

0.45331 

.800 

18 

0.60275 

0.52695 

0.40069 

0.21591 

0.47832 

.800 

19 

0.62261 

0.55500 

0.42086 

0.22006 

0.50334 

.800 

20 

0.64148 

0.58306 

0.44102 

0.22402 

0.52836 

H  — — 

means  value  is  higher  than  that 
independence . 

given  by 

conservative 

assumption  of 
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Table  la  Continued. 


2=1.96 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariate 

=  P.  . 
il  iJ 

ij  Jj  ’eT’ 

Equicorrelated 

Case  P,  =  P 
k£ 

P..,  =  (P..)2 
ij  ij 

ij.jj'eT' 

all  k  ,4 

.900 

3 

0.14261 

0.09070 

0.08364 

0.08364 

0.08764 

.900 

4 

0.18548 

0.11105 

0.09693 

0.09359 

0.10493 

.900 

5 

0.22620 

0.13141 

0.11023 

0.10157 

0.12223 

.900 

6 

0.26489 

0.15176 

0.12352 

0.10825 

0.13952 

.900 

7 

0.30164 

0.17211 

0.13681 

0.11399 

0.15681 

.900 

8 

0.33656 

0.19246 

0.15011 

0.11904 

0.17410 

.900 

9 

0.36973 

0.21282 

0.16340 

0.12354 

0.19140 

.900 

10 

0.40124 

0.23317 

0.17669 

0.12761 

0.20869 

.900 

11 

0.43117 

0.25352 

0.18999 

0.13132 

0.22598 

.900 

12 

0.45961 

0.27387 

0.20328 

0.13473 

0.24328 

.900 

13 

0.48663 

0.29423 

0.21657 

0.13788 

0.26057 

.900 

14 

0.51229 

0.31458 

0.22987 

0.14082 

0.27786 

.900 

15 

0.53668 

0.33493 

0.24316 

0.14357 

0.29515 

.900 

16 

0.55984 

0.35528 

0.25645 

0.14615 

0.31245 

.900 

17 

0.58185 

0.37564 

0.26975 

0.14859 

0.32974 

.900 

18 

0.60275 

0.39599 

0.28304 

0.15090 

0.34703 

.900 

19 

0.62261 

0.41634 

0.29633 

0.15309 

0.36433 

.900 

20 

0.64148 

0.43669 

0.30963 

0.15517 

0.38162 

.990 

3 

0.14261 

0,06315 

0.06018 

0.06018 

* 

.990 

4 

0.18548 

0.06973 

0.06378 

0.06264 

* 

.990 

5 

0.22620 

0.07631 

0.06739 

0.06449 

★ 

.990 

6 

0.26489 

0.08289 

0.07099 

0.06597 

* 

.990 

7 

0.30164 

0.08947 

0.07460 

0.06720 

* 

.990 

8 

0.33656 

0.09605 

0.07820 

0.06825 

* 

.990 

9 

0.36973 

0.10263 

0.08181 

0.06917 

* 

.990 

10 

0.40123 

0.10920 

0.08541 

0.06998 

* 

.990 

11 

0.43117 

0.11578 

0.08902 

0.07071 

* 

.990 

12 

0.45961 

0.12236 

0.09262 

0.07136 

* 

.990 

13 

0.48663 

0.12894 

0.09623 

0.07196 

k 

.990 

14 

0.51229 

0.13552 

0.09982 

0.07252 

k 

.990 

15 

0.53668 

0.14210 

0.10344 

0.07303 

* 

.990 

16 

0.55984 

0.14868 

0.10704 

0.07350 

k 

.990 

17 

0.58185 

0.15526 

0.11065 

0.07395 

k 

.990 

18 

0.60275 

0.16183 

0.11425 

0.07436 

k 

.990 

19 

0.62261 

0.16841 

0.11786 

0.07475 

k 

.990 

20 

0.64148 

0.17499 

0.12146 

0.07513 

k 

m 


i 


w 

I 


II 


means  value  is  missing  due  to  inability  of  LMSL  subroutine  DJ'lLIN 
to  integrate. 
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Table  lb.  Upper  bounds  for  Pr { (  X .  )£ (-z ,z) } 

i=l  1 


Mfi* 

1 


2=2.25 

RHO  Number  Independence  Bivariate  Trivariate  Exact  for 


Trivariate 


P..,  =  P..  Equicorrelated  P . . ,  =  (P..)' 

«  Case  P  s  P  «  U 

ijjJJ  eT  ..f 

all  k,2.  13,33 


0.07137 


.100 

3 

0.07157 

0.07195 

0.07128 

0.07128 

.100 

4 

0.09427 

— 

— 

0.09455 

.100 

5 

0.11641 

— 

— 

0.11625 

.100 

6 

0.13801 

— 

— 

0.13726 

.100 

7 

0.15909 

— 

— 

0.15789 

.100 

8 

0.17965 

— 

— 

0.17786 

.100 

9 

0.19971 

— 

— 

0.19730 

.100 

10 

0.21927 

— 

— 

0.21622 

.100 

11 

0.23836 

— 

— 

0.23464 

.100 

12 

0.25698 

— 

— 

0.25256 

.100 

13 

0.27515 

— 

— 

0.27001 

.100 

14 

0.29287 

— 

— 

0.28701 

.100 

15 

0.31016 

— 

— 

0.30357 

.100 

16 

0.32702 

— 

— — 

0.31969 

.100 

17 

0.34348 

— 

— 

0.33540 

.100 

18 

0.35953 

— 

— 

0.35070 

.100 

19 

0.37519 

— 

0.36561 

.100 

20 

0.39046 

— 

— 

0.38010 

.200 

3 

0.07157 

0.07132 

0.07039 

0.07039 

.200 

4 

0.09427 

— 

0.09290 

0.09318 

.200 

5 

0.11641 

— 

0.11540 

0.11380 

.200 

6 

0.13801 

— 

0.13791 

0.13370 

.200 

7 

0.15909 

— 

— 

0.15242 

.200 

8 

0.17965 

— 

— 

0.17151 

.200 

9 

0.19971 

— 

— 

0.18949 

.200 

10 

0.21927 

— 

— 

0.20691 

.200 

11 

0.23836 

— 

— 

0.22377 

.200 

12 

0.25698 

— 

— 

0.24012 

.200 

13 

0.27515 

— 

— 

0.25598 

.200 

14 

0.29287 

— 

— 

0.27136 

.200 

15 

0.31016 

— 

— 

0.28630 

.200 

16 

0.32702 

— 

— 

0 . 30080 

.200 

17 

0.34348 

— 

— 

0.31490 

.200 

18 

0.35953 

— 

— 

0.32861 

.200 

19 

0.37519 

— 

— 

0.34194 

.200 

20 

0.39046 

— 

— 

0.35491 

*„ _ 

"  means 

value  is  higher 

than  that 

given  by 

conservative 

0.07075 

0.09362 
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independence. 
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Table  lb  Continued. 


2=2.25 

RHO  Number  Independence  Bivariate 


.300 

3 

0.07157 

0.07025 

.300 

4 

0.09427 

0.09315 

.300 

5 

0.11641 

0.11606 

.300 

6 

0.13801 

— 

.300 

7 

0.15909 

— 

.300 

8 

0.17965 

— 

.300 

9 

0.19971 

— 

.300 

10 

0.21927 

— 

.300 

11 

0.23836 

— 

.300 

12 

0.25698 

— 

.300 

13 

0.25715 

— 

.300 

14 

0.29287 

— 

.300 

15 

0.31016 

— 

.300 

16 

0.32702 

— 

.300 

17 

0.34348 

— 

.300 

18 

0.35953 

— 

.300 

19 

0.37519 

0.43668 

.300 

20 

0.39046 

0.45958 

.400 

3 

0.07157 

0.06363 

.400 

4 

0.09427 

0.09079 

.400 

5 

0.11641 

0.11291 

.400 

6 

0.13801 

0.13502 

.400 

7 

0.15909 

0.15714 

.400 

3 

0.17965 

0.17925 

.400 

9 

0.19971 

— 

.400 

10 

0,21927 

— 

.400 

11 

0.23836 

— 

.400 

12 

0.25698 

— 

.400 

13 

0.27515 

— 

.400 

14 

0.29287 

— 

.400 

15 

0.31016 

— 

.400 

16 

0.32702 

— 

.400 

17 

0.34348 

— 

.400 

18 

0.35953 

— 

.400 

19 

0.37519 

— 

.400 

20 

0.39046 

— 

Trivariate 

Exact  for 

Trivariate 

P, . ,  =  P.  . 
ij  iJ 

ij.jj’eT' 

Equicorrelated 

Case  P.  .  s  P 
k£ 

V  '  (pi3 

ij ,jj 'eT1 

all  k,  £ 

0.06892 

0.06892 

0.06968 

0.09049 

0.08572 

0.09200 

0.11205 

0.10886 

0.11433 

0.13362 

0.12723 

0.13665 

0.15519 

0.14478 

0.15898 

0.17676 

0.16158 

— 

0.19833 

0.17770 

— 

— 

0.19318 

— 

— 

0.20808 

— 

— 

0.22243 

— 

— 

0.23628 

— 

— 

0.24966 

— 

— 

0.26259 

— 

— - 

0.27511 

— 

— 

0.28724 

— 

— 

0.29900 

— 

0.41401 

0.31041 

— 

0.43557 

0.32150 

— 

0.06684 

0.06684 

0.06306 

0.08711 

0.08572 

0.08955 

0.1073S 

0.10317 

0.11105 

0.12766 

0.11956 

0.13254 

0.14793 

0.13503 

0.15404 

0.16820 

0.14968 

0.17553 

0.13843 

0.16359 

0.19702 

0.20875 

0.17686 

0.21852 

0.22902 

0.13953 

— 

0.24930 

0.20167 

— 

0.26957 

0.21332 

— 

0.28984 

0.22452 

— 

0.31011 

0.23530 

— 

— 

0.24571 

— 

— 

0.25575 

— 

— 

0.26547 

— 

— 

0.27487 

— 

— 

0.28399 

— 

means  value  is  higher  than  that  given  by  conservative  assumption  of 
independence . 


Table  lb  Continued. 


I 


m 

2*2.25 

RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariate 

M 

P. *  P.  . 
il  1J 

ij.jl'eT’ 

Equicorrelated 

P<.,  = 

r» 

Case  P.  „  a  P 
k2. 

all  k,jt 

ij  ij 

ijJj’eT' 

as 

i 

Vi. 

.500 

3 

0.07157 

0.06651 

0.06410 

0.06410 

0.06579 

w 

iVtV 

.500 

4 

0.09427 

0.08754 

0.08272 

0.08090 

0.08610 

.500 

5 

0.11641 

0.10857 

0.10134 

0.09625 

0.10641 

if 

.500 

6 

0.13801 

0.12960 

0.11996 

0.11042 

0.12672 

.500 

7 

0.15909 

0.15063 

0.13858 

0.12361 

0.14703 

;Jtiv 

.500 

8 

0.17965 

0.17166 

0.15720 

0.13595 

0.16734 

M 

.500 

9 

0.19971 

0.19269 

0.17582 

0.14757 

0.18765 

k 

I 

i 

.500 

10 

0.21927 

0.21372 

0.19444 

0.15855 

0.20796 

.500 

11 

0.23836 

0.23475 

0.21306 

0.16896 

0.22827 

.500 

12 

0.25698 

0.25578 

0.23168 

0.17887 

0.24858 

.500 

13 

0.27515 

— 

0.25030 

0.18833 

0.26889 

.500 

14 

0.29287 

— 

0.26892 

0.19739 

0.28921 

.500 

15 

0.31016 

— 

0.28754 

0.20607 

0.30952 

.500 

16 

0.32702 

— 

0.30616 

0.21441 

— 

.500 

17 

0.34343 

— 

0.32478 

0.22433 

— 

i 

.500 

18 

0.35953 

— 

0.34340 

0.23017 

— 

.500 

19 

0.37519 

— 

0.36202 

0.23765 

— 

P 

.500 

20 

0.39046 

— 

0.38064 

0.24487 

— 

.600 

3 

0.07157 

0.06359 

0.06061 

0.06061 

0.06271 

.600 

4 

0.09427 

0.08316 

0.07720 

0.07513 

0.08140 

S*l& 

.600 

5 

0.11641 

0.10273 

0.09379 

0.08811 

0.10009 

i 

.600 

6 

0.13801 

0.12230 

0.11038 

0.09934 

0.11878 

.600 

7 

0.15909 

0.14183 

0.12697 

0.11068 

0.13747 

m 

.600 

8 

0.17965 

0.16145 

0.14356 

0.12067 

0.15616 

.600 

9 

0.19971 

0.18102 

0.16015 

0.12998 

0.17485 

_u 

.600 

10 

0.21927 

0.20059 

0.17674 

0.13870 

0.19355 

ij® 

.600 

11 

0.23G36 

0.22016 

0.19333 

0.14692 

0.21224 

.600 

12 

0.25698 

0.23973 

0.20992 

0.15469 

0.23093 

I 

.600 

13 

0.27515 

0.25930 

0.22651 

0.16207 

0.24962 

.600 

14 

0.29287 

0.27887 

0.24309 

0.16909 

0.26831 

.600 

15 

0.31016 

0.29844 

0.25968 

0.17580 

0.28700 

n 

.600 

16 

0.32702 

0.31802 

0.27627 

0.18221 

0.30569 

SSt 

.600 

17 

0.34348 

0.33759 

0.29286 

0.18837 

0.32438 

.600 

18 

0.35953 

0.35716 

0.30945 

0.19428 

0.34307 

.600 

19 

0.37519 

— 

0.32604 

0.19997 

0.36176  1 

m 

.600 

20 

0.39046 

— 

0.34263 

0.20547 

0.38045 

" - "  means  value  is  higher  than  that  given  by  conservative  assumption  of 

independence. 
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Table  lb  Continued. 


z=2. 25 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariati 

P  =  p 

i  j  ’  ij 

ij.jj'eT’ 

Equicorrelated 
Case  Pu  .  P 

V  *  (p 

iJ.Jj'eT' 

all  k,£ 

.700 

3 

0.07157 

0.05968 

0.05620 

0.05620 

0.05858 

.70  0 

4 

0.09427 

0.07729 

0.07033 

0.06814 

0.07511 

.700 

5 

0.11641 

0.09491 

0.08447 

0.07854 

0.09163 

.700 

6 

0.13801 

0.11252 

0.09860 

0.08780 

0.10815 

.700 

7 

0.15909 

0.13013 

0.11274 

0.09617 

0.12467 

.700 

8 

0.17965 

0.14775 

0.12687 

0.10381 

0.14119 

.700 

9 

0.19971 

0.16536 

0.14101 

0.11086 

0.15771 

.700 

10 

0.21927 

0.18298 

0.15514 

0.11741 

0.17424 

.700 

11 

0.23836 

0.20059 

0.16927 

0.12353 

0.19076 

.700 

12 

0.25698 

0.21820 

0.18341 

0.12928 

0.20728 

.700 

13 

0.27515 

0.23582 

0.19754 

0.13471 

0.22380 

.700 

14 

0.29287 

0.25343 

0.21168 

0.13985 

0.24032 

.700 

15 

0.31016 

0.27105 

0.22581 

0.14473 

0.25684 

.700 

16 

0.32702 

0.28866 

0.23995 

0.14934 

0.27337 

.700 

17 

0.34348 

0.30628 

0.25408 

0.15383 

0.28989 

.700 

18 

0.35953 

0.32389 

0.26822 

0.15809 

0.30641 

.700 

19 

0.37519 

0.34150 

0.28235 

0.16218 

0.32293 

.700 

20 

0.39046 

0.35912 

0.29649 

0.16611 

0.33945 

.800 

3 

0.07157 

0.05430 

0.05053 

0.05053 

0.05299 

.800 

4 

0.09427 

0.06923 

0.06169 

0.05956 

0.06659 

.800 

5 

0.11641 

9,^8416 

0.07284 

0.06714 

0.08020 

.800 

6 

0.13801 

0.09909 

0.08399 

0.07384 

0.09381 

.800 

7 

0.15909 

0.11401 

0.09515 

0.07974 

0.10742 

.800 

8 

0.17965 

0.12895 

0.10630 

0.08505 

0.12103 

.800 

9 

0.19971 

0.14387 

0.11746 

0.08990 

0.13464 

.800 

10 

0.21927 

0.15880 

0.12861 

0.09439 

0.14825 

.800 

11 

0.23836 

0.17372 

0.13977 

0.09848 

0.16186 

.800 

12 

0.25698 

0.18865 

0.15092 

0.10233 

0.17547 

.800 

13 

0.27515 

0.20358 

0.16208 

0.10594 

0.18908 

.800 

14 

0.29287 

0.21851 

0.17323 

0.10934 

0.20268 

.800 

15 

0.31016 

0.23343 

0.18439 

0.11255 

0.21629 

.800 

16 

0.32702 

0.24836 

0.19554 

0.11560 

0.22990 

.800 

17 

0.34348 

0.26329 

0.20670 

0.11850 

0.24351 

.800 

18 

0.35953 

0.27821 

0.21785 

0.12126 

0.25712 

.800 

19 

0.37519 

0.29314 

0.22900 

0.12391 

0.27073 

.800 

20 

0.39046 

0.30807 

0.24016 

0.12644 

0.28434 

ij 
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Table  lb  Continued. 


z-2 . 25 


RHO  Number  Independence  Bivariate  Trivariate  Exact  for 


Trivariate 


P..,  *  P.. 

Equicorrelated 

P..,  =  (P 

i3  13 

Case  P,  =  P 

ij.jj’eT* 

k£  _ 

ij.jj’eT’ 

all  k,£ 

.900 

3 

0.07157 

0.04633 

0.04277 

0.04277 

0.04485 

.900 

4 

0.09427 

0.05727 

0.05015 

0.04840 

0.05437 

.900 

5 

0.11641 

0.06821 

0.05753 

0.05299 

0.06375 

.900 

6 

0.13801 

0.07915 

0.06491 

0.05686 

0.07232 

.900 

7 

0.15909 

0.09010 

0.07229 

0.06022 

0.08278 

.900 

8 

0.17965 

0.10104 

0.07967 

0.06319 

0.09231 

.900 

9 

0.19971 

0.11198 

0.08705 

0.06586 

0.10752 

.900 

10 

0.21927 

0.12292 

0.09443 

0.06820 

0.11019 

.900 

11 

0.23836 

0.13386 

0.10181 

0.07051 

0.12132 

.900 

12 

0.25698 

0.14480 

0.10910 

0.07256 

0.13030 

.900 

13 

0.27515 

0.15574 

0.11657 

0.07447 

0.13979 

.900 

14 

0.29287 

0.16668 

0.12395 

0.07625 

0.14882 

.900 

15 

0.31016 

0.17762 

0.13133 

0.07792 

0.15903 

.900 

16 

0.32702 

0.18856 

0.13871 

0.07950 

0.16805 

.900 

17 

0.34348 

0.19951 

0.14609 

0.08099 

0.17794 

.900 

18 

0.35953 

0.21045 

0.15347 

0.08241 

0.18621 

.900 

19 

0.37519 

0.22139 

0.16085 

0.08375 

0.19728 

.900 

20 

0.39046 

0.23233 

0.16823 

0.08504 

0.20571 

.990 

3 

0.07157 

0.03159 

0.03000 

0.03000 

* 

.990 

4 

0.09427 

0.03516 

0.03198 

0.03136 

* 

.990 

5 

0.11641 

0.03873 

0.03396 

0.03234 

* 

.990 

6 

0.13801 

0.04230 

0.03594 

0.03322 

* 

.990 

7 

0.15909 

0.04586 

0.03792 

0.03391 

* 

.990 

8 

0.17965 

0.04943 

0.03990 

0.03450 

* 

.990 

9 

0.19971 

0.05300 

0.04188 

0.03511 

* 

.990 

10 

0.21927 

0.05657 

0.04387 

0.03547 

* 

.990 

11 

0.23836 

0.06014 

0.04585 

0.03588 

* 

.990 

12 

0.25698 

0.06371 

0.04783 

0.03625 

* 

.990 

13 

0.27515 

0.06728 

0.04981 

0.03659 

* 

.990 

14 

0.29287 

0.07085 

0.05179 

0.03691 

* 

.990 

15 

0.31016 

0.07442 

0.05377 

0.03719 

* 

.990 

16 

0.32702 

0.07799 

0.05575 

0.03746 

* 

.990 

17 

0.34348 

0.08156 

0.05773 

0.03771 

* 

.990 

18 

0.35935 

0.08513 

0.05971 

0.03795 

* 

.990 

19 

0.37519 

0.08870 

0.06169 

0.03817 

* 

.990 

20 

0.39046 

0.09227 

0.06368 

0.03834 

* 

"*"  means  value  is  missing  due  to  the  inability  of  DMLIN  to  integrate. 
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Table  lc.  Upper  bounds  for  Pr{ (  X. )^(-z,z) } 

i*l  X 

z=2.50 

RHO  Number  Independence  Bivariate  Trivariate  Exact  for  Trivariate 


.100 

3 

0.03680 

— 

.100 

4 

0.04876 

— 

.100 

5 

0.06057 

— 

.100 

6 

0.07224 

— 

.100 

7 

0.08376 

— 

.100 

8 

0.09514 

— 

.100 

9 

0.10638 

— 

.100 

10 

0.11748 

— 

.100 

11 

0.12844 

— 

.100 

12 

0.13926 

— 

.100 

13 

0.14995 

— 

.100 

14 

0.16051 

— 

.100 

15 

0.17093 

— 

.100 

16 

0.18123 

— 

.100 

17 

0.19140 

— 

.100 

18 

0.20144 

— 

.100 

19 

0.21136 

— 

.100 

20 

0.22115 

- — 

.200 

3 

0.03680 

0.03663 

.200 

4 

0.04876 

0.04874 

.200 

5 

0.06057 

— 

.200 

6 

0.07224 

— 

.200 

7 

0.08376 

— 

.200 

8 

0.09514 

— 

.200 

9 

0.10638 

— 

.200 

10 

0.11748 

— 

.200 

11 

0.12844 

— 

.200 

12 

0.13926 

— 

.200 

13 

0.14995 

— 

.200 

14 

0.16051 

— 

.200 

15 

0.17093 

— 

.200 

16 

0.18123 

— 

.200 

17 

0.19140 

— 

.200 

18 

0.20144 

— 

.200 

19 

0.21136 

— 

.200 

20 

0.22115 

— 

" - "means  value  is  higher  than  that 

independence. 


P...  =  P..  Equicorrelated  P . . ,  =  (P..) 

1J  1J  Case  P  =  P  ij  i2 

ij.jj'eT'  °aSe  -  ij.jj'eT' 

all  k.,2. 


0.03669 

0.03669 

0.03672 

0.04873 

0.04854 

— 

— 

0.06220 

— 

— 

0.07170 

— 

— 

0.08303 

— 

— 

0.09418 

— 

— 

0.10517 

— 

— 

0.11599 

— 

— 

0.12665 

— 

— 

0.13716 

— 

— 

0.14752 

— 

— 

0.15772 

— 

— 

0.16778 

— 

— 

0.17769 

— 

— 

0.18746 

— 

— 

0.19710 

— 

— 

0.20659 

— 

— 

0.21595 

— 

0.03634 

0.03634 

0.03648 

0.04815 

0.04786 

0.04844 

0.05996 

0.05913 

0.06039 

0.07177 

0.07013 

— 

0.08359 

0.09090 

— 

— 

0.09143 

— 

— 

0.10516 

— 

— 

0.11599 

— 

— 

0.12174 

— 

— 

0.13144 

— 

— 

0.14095 

— 

- — 

0.15028 

— 

— 

0.15943 

— 

0.16288 

— 

— 

0.17724 

— 

— 

0.18590 

— 

— 

0.19441 

— 

— 

0.22781 

— 

given  by  conservative  assumption  of 


Table  lc  Continued 


z*2 . 50 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

P..,  =  P. . 
ij  i  3 

ij.jj’eT' 

Exact  for 
Equicorrelated 
Case  Pu  .  P 

all  k,£ 

Trivariate 

P. =  (P. 

x. 

ij.jj'eT' 

.300 

3 

0.03680 

0.03620 

0.03573 

0.03573 

0.03604 

.300 

4 

0.04876 

0.04809 

0.04715 

0.04673 

0.04777 

.300 

5 

0.06057 

0.05999 

0.05857 

0.05733 

0.05950 

.300 

6 

0.07224 

0.07188 

0.06999 

0.06758 

0.07132 

.300 

7 

0.08376 

— 

0.08142 

0.07750 

0.08296 

.300 

8 

0.09514 

— 

0.09284 

0.08712 

0.09470 

.300 

9 

0.10638 

— 

0.10426 

0.09646 

— 

.300 

10 

0.11748 

— 

0.11568 

0.10553 

— 

.300 

11 

0.12844 

— 

0.12710 

0.11436 

— 

.300 

12 

0.13926 

— 

0.13852 

0.12295 

— 

.300 

13 

0.14995 

— 

0.14994 

0.13133 

— 

.300 

14 

0.16051 

— 

— 

0.13950 

— 

.300 

15 

0.17093 

— 

— 

0.14747 

— 

.300 

16 

0.18123 

— 

— 

0.15526 

— 

.300 

17 

0.19140 

— 

— 

0.16288 

— 

.300 

18 

0.20144 

— 

— 

0.17032 

— 

.300 

19 

0.21136 

— 

— 

0.17761 

— 

.300 

20 

0.22115 

— 

— 

0.18475 

— 

.400 

3 

0.03680 

0.03554 

0.03483 

0.03485 

0.03535 

.400 

4 

0.04876 

0.04710 

0.04568 

0.04509 

0.04672 

.400 

5 

0.06057 

0.05866 

0.05653 

0.05481 

0.05809 

.400 

6 

0.07224 

0.07022 

0.06738 

0.06408 

0.06946 

.400 

7 

0.08376 

0.08178 

0.07824 

0.07295 

0.08083 

.400 

8 

0.09514 

0.09334 

0.08909 

0.08145 

0.09220 

.400 

9 

0.10638 

0.10491 

0.09994 

0.08962 

0.10357 

.400 

10 

0.11748 

0.11647 

0.11079 

0.09749 

0.11494 

.400 

11 

0.12844 

0.12803 

0.12164 

0.10509 

0.12631 

.400 

12 

0.13926 

— 

0.13249 

0.11244 

0.13768 

.400 

13 

0.14995 

— 

0.14334 

0.11956 

0.14905 

.400 

14 

0.16051 

— 

0.15419 

0.12645 

0.16043 

.400 

15 

0.17093 

— 

0.16504 

0.13315 

— 

.400 

16 

0.18123 

— 

0.17590 

0.13966 

— 

.400 

17 

0.19140 

— 

0.18675 

0.14599 

— 

.400 

18 

0.20144 

— 

0.19760 

0.15216 

— 

.400 

19 

0.21136 

— 

0.20845 

0.15812 

— 

.400 

20 

0.22115 

— 

0.21930 

0.16403 

— 

- "  means  value  is  higher  than  that  given  by  conservative  assumption  of 

independence. 


Table  lc  Continued. 
z»2.50 

RHO  Number  Independence  Bivariate  Trivariate 

Pij*  =  Pij 
ijjj'tr 


Exact  for 

Equicorrelated 

Case  P,  .  s  P 
k£ 

all  k.S. 


Trivariate 

V  '  (V 

ij.jj’cT’ 


.500 

3 

0.03680 

0.03458 

.500 

4 

0.04876 

0.04566 

.500 

5 

0.06057 

0.05674 

.500 

6 

0.07224 

0.06782 

.500 

7 

0.08376 

0.07890 

.500 

8 

0.09514 

0.08998 

.500 

9 

0.10638 

0.10106 

.500 

10 

0.11748 

0.11214 

.500 

11 

3.12844 

0.12322 

.500 

12 

0.13926 

0.13431 

.500 

13 

0.14995 

0.14539 

.500 

14 

0.16051 

0.15647 

.500 

15 

0.17093 

0.16755 

.500 

16 

0.18123 

0.17863 

.500 

17 

0.19140 

0.18971 

.500 

18 

0.20144 

0.20079 

.500 

19 

0.21136 

— 

.500 

20 

0.22115 

— 

.600 

3 

0.03680 

0.03322 

.600 

4 

0.04876 

0.04362 

.600 

5 

0.06057 

0.05402 

.600 

6 

0.07224 

0.06443 

.600 

7 

0.08376 

0.07483 

.600 

8 

0.09514 

0.08523 

.600 

9 

0.10638 

0.09563 

.600 

10 

0.11748 

0.10603 

.600 

11 

0.12844 

0.11643 

.600 

12 

0.13926 

0.12683 

.600 

13 

0.14995 

0.13723 

.600 

14 

0.16051 

0.14764 

.600 

15 

0.17093 

0.15804 

.600 

16 

0.18123 

0.16844 

.600 

17 

0.19140 

0.17884 

.600 

18 

0.20144 

0.18924 

.600 

19 

0.21136 

0.19964 

.600 

20 

0.22115 

0.21004 

" - "  means  value  is  higher  than 

independence. 


0.03358 

0.03358 

0.03433 

0.04367 

0.04288 

0.04517 

0.05375 

0.05154 

0.05600 

0.06383 

0.05965 

0.06684 

0.07391 

0.06731 

0.07767 

0.08400 

0.07456 

0.08851 

0.09408 

0.08146 

0.09934 

0.10416 

0.08805 

0.11018 

0.11425 

0.09436 

0.12101 

0.12433 

0.10042 

0.13184 

0.13441 

0.10624 

0.14268 

0.14450 

0.11856 

0.15351 

0.15458 

0.12028 

0.16435 

0.16466 

0.12252 

0.17518 

0.17475 

0.12760 

0.18602 

0.18483 

0.13252 

0.19685 

0.19491 

0.13706 

0.20769 

0.20499 

0.14195 

0.21852 

0.03191 

0.03191 

0.03289 

0.04100 

0.04004 

0.04296 

0.05008 

0.04745 

0.05302 

0.05917 

0.05427 

0.06309 

0.06826 

0.06060 

0.07316 

0.07734 

0.06653 

0.08323 

0.08643 

0.07211 

0.09330 

0.09552 

0.07739 

0.10337 

0.10461 

0.08240 

0.11343 

0.11369 

0.08718 

0.12350 

0.12278 

0.09174 

0.13357 

0.13187 

0.09611 

0.14364 

0.14095 

0.10031 

0.15371 

0.15004 

0.10435 

0.16377 

0.15913 

0.10825 

0.17384 

0.16822 

0.11201 

0.18391 

0.17730 

0.11565 

0.19398 

0.18639 

0.11918 

0.20405 

given  by 

conservative  , 

assumption  of 

/«  .'N 


Table  lc  Continued. 


z-2.50 

RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariat 

.700 

3 

0.03680 

0.03132 

P,, ,  =  P.  . 
ij  13 

ij.jj’eT' 

0.02969 

Equicorrelated 
Case  P,  .  s  P 

all  k ,  2. 

0.02969 

Pij •  =  (P 
ij.jj’eT* 

0.03086 

.700 

4 

0.04876 

0.04076 

0.03752 

0.03645 

0.03986 

.700 

5 

0.06057 

0.05021 

0.04534 

0.04244 

0.04885 

.700 

6 

0.07224 

0.05966 

0.05317 

0.04785 

0.05785 

.700 

7 

0.08376 

0.06911 

0.06099 

0.05279 

0.06684 

.700 

8 

0.09514 

0.07855 

0.06882 

0.05735 

0.07584 

.700 

9 

0.10638 

0.08800 

0.07664 

0.06159 

0.08483 

.700 

10 

0.11748 

0.09745 

0.08447 

0.06557 

0.09383 

.700 

11 

0.12844 

0.10690 

0.09229 

0.06930 

0.10283 

.700 

12 

0.13926 

0.11635 

0.10011 

0.07284 

0.11182 

.700 

13 

0.14995 

0.12579 

0.10794 

0.07620 

0.12082 

.700 

14 

0.16051 

0.13524 

0.11576 

0.07939 

0.12981 

.700 

15 

0.17093 

0.14469 

0.12359 

0.08244 

0.13881 

.700 

16 

0.18123 

0.15415 

0.13141 

0.08536 

0.14780 

.700 

17 

0.19140 

0.16359 

0.13924 

0.08817 

0.15680 

.700 

18 

0.20144 

0.17303 

0.14706 

0.09087 

0.16579 

.700 

19 

0.21136 

0.18248 

0.15489 

0.09346 

0.17479 

.700 

20 

0.22115 

0.19193 

0.16271 

0.09597 

0.18378 

.800 

3 

0.03680 

0.02859 

0.02673 

0.02673 

0.02799 

.800 

4 

0.04876 

0.03667 

0.03296 

0.03188 

0.03548 

.800 

5 

0.06057 

0.04475 

0.03919 

0.03629 

0.04297 

.800 

6 

0.07224 

0.05284 

0.04542 

0.04018 

0.05047 

.800 

7 

0.08376 

0.06092 

0.05165 

0.04367 

0.05796 

.800 

8 

0.09514 

0.06900 

0.05788 

0.04684 

0.06545 

.800 

9 

0.10638 

0.07709 

0.06411 

0.04988 

0.07294 

.800 

10 

0.11748 

0.08517 

0.07034 

0.05244 

0.08043 

.800 

11 

0.12844 

0.09325 

0.07657 

0.05500 

0.08792 

.800 

12 

0.13926 

0.10134 

0.08280 

0.05731 

0.09541 

.800 

13 

0.14995 

0.10942 

0.08903 

0.05952 

0.10290 

.800 

14 

0.16051 

0.11750 

0.09526 

0.06162 

0.11039 

.800 

15 

0.17093 

0.12559 

0.10149 

0.06361 

0.11788 

.800 

16 

0.18123 

0.13367 

0.10772 

0.06550 

0.12537 

.800 

17 

0.19140 

0.14175 

0.11395 

0.06731 

0.13286 

.800 

18 

0.20144 

0.14984 

0.12018 

0.06904 

0.14035 

.800 

19 

0.21136 

0.15792 

0.12641 

0.07070 

0.14784 

.800 

20 

0.22115 

0.16600 

0.13264 

0.07224 

0.15534 
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.900 

3 

0.03680 

0.02439 

0.02255 

0.02255 

0.02368 

.900 

4 

0.04876 

0.03037 

0.02669 

0.02576 

0.02895 

.900 

5 

0.06057 

0.03635 

0.03083 

0.02840 

0.03422 

.900 

6 

0.07224 

0.04233 

0.03498 

0.03066 

0.03950 

.900 

7 

0.08376 

0.04832 

0.03912 

0.03263 

0.04477 

.900 

8 

0.09514 

0.05430 

0.04326 

0.03439 

0.05004 

.900 

9 

0.10638 

0.06028 

0.04741 

0.03597 

0.05531 

.900 

10 

0.11748 

0.06626 

0.05155 

0.03742 

0.06059 

.900 

11 

0.12844 

0.97225 

0.05569 

0.03875 

0.06586 

.900 

12 

0.13926 

C.u. 323 

0.05984 

0.03998 

0.07113 

.900 

13 

0.14995 

0.08421 

0.06398 

0.04113 

0.07641 

.900 

14 

0.16051 

0.09020 

0.06813 

0.04211 

0.08168 

.900 

15 

0.17093 

0.09618 

0.07227 

0.04322 

0.08695 

.900 

16 

0.18123 

0.10216 

0.07641 

0.04417 

0.09223 

.900 

17 

0.19140 

0.10814 

0.08056 

0.04508 

0.09750 

.900 

18 

0.20144 

0.11413 

0.08470 

0.04594 

0.10277 

.900 

19 

0.21136 

0.12011 

0.08884 

0.04677 

0.10805 

.900 

20 

0.22115 

0.12609 

0.09299 

0.04755 

0.11332 

.990 

3 

0.03680 

0.01636 

0.01549 

0.01549 

* 

.990 

4 

0.04876 

0.01833 

0.01660 

0.01626 

* 

.990 

5 

0.06057 

0.02030 

0.01771 

0.01684 

* 

.990 

6 

0.07224 

0.02227 

0.01881 

0.01730 

* 

.990 

7 

0.08376 

0.02423 

0.01991 

0.01770 

* 

.990 

8 

0.09514 

0.02620 

0.02102 

0.01803 

* 

.990 

9 

0.10638 

0.02817 

0.02212 

0.01832 

* 

.990 

10 

0.11748 

0.03014 

0.02323 

0.01858 

* 

.990 

1.1 

0.12844 

0.03211 

0.02433 

0.01881 

* 

.990 

12 

0.13926 

0.03408 

0.02544 

0.01903 

"k 

.990 

13 

0.14995 

0.03605 

0.02654 

0.01922 

* 

.990 

14 

0.16051 

0.03802 

0.02765 

0.01940 

* 

.990 

15 

0.17093 

0.03999 

0.02875 

0.01957 

* 

.990 

16 

0.18123 

0.04196 

0.02986 

0.01972 

* 

.990 

17 

0.19140 

0.04393 

0.03096 

0.01987 

* 

.990 

18 

0.20144 

0.04590 

0.03207 

0.02000 

* 

.990 

19 

0.21136 

0.04787 

0.03317 

0.02013 

* 

.990 

20 

0.22115 

0.04983 

0.03428 

0.02025 

* 

* 

means 

value  is  missing 

due  to  inability  of 

DMLIN  to  integrate. 

Table  Id.  Upper  Bounds  for  Pr{(  (J  X.)  i  (-z,z)} 

i-1  1 


z=2.75 

RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariate 

.100 

3 

0.01777 

P.  =  P. . 

iJ  13 

ij.jj’eT’ 

0.01774 

Equicorrelated 

Case  P,  „  s  P 
k£ 

all  k,£ 

0.01774 

P. .  =  (P.  . 

13  13 

ij.jj’eT' 

0.01775 

.100 

4 

0.02363 

— 

0.02361 

0.02359 

0.02363 

.100 

5 

0.02944 

— 

— 

0.02945 

_ 

.100 

6 

0.03523 

— 

— 

0.03517 

_ 

.100 

V 

0.04098 

— 

— 

0.04084 

_ 

.100 

8 

0.04669 

— 

— 

0.04647 

_ 

.100 

9 

0.05237 

— 

— 

0.05205 

_ 

.100 

10 

0.05802 

— 

— 

0.05759 

_ 

.100 

11 

0.06364 

— 

— 

0.06309 

_ 

.10C 

12 

0.06922 

— 

— 

0.06855 

_ 

.100 

13 

0.07476 

— 

— 

0.07396 

- — 

.100 

14 

0.08028 

— 

— 

0.07933 

— 

.100 

15 

0.08576 

— 

— 

0.08467 

— 

.100 

16 

0.09121 

— 

— 

0.08995 

— 

.100 

17 

0.09662 

— 

— 

0.09521 

— 

.100 

18 

0.10201 

— 

— 

0.10042 

— 

.100 

19 

0.10736 

— 

— 

0.10559 

— 

.100 

20 

0.11268 

— 

— 

0.11073 

— 

.200 

3 

0.01777 

0.01770 

0.01762 

0.01762 

0.01767 

.200 

4 

0.02363 

0.02357 

0.02340 

0.02341 

0.02350 

.200 

5 

0.02944 

0.02944 

0.02919 

0.02945 

0 . 02934 

.200 

6 

0.03523 

— 

0.03498 

0.03455 

0.03518 

.200 

7 

0.04098 

— 

0.04077 

0.04002 

— 

.200 

8 

0.04669 

— 

0.04655 

0.04541 

— 

.200 

9 

0.05237 

— 

0.05234 

0.05073 

— 

.200 

10 

0.05802 

— 

— 

0.05599 

— 

.200 

11 

0.06364 

— 

— 

0.06118 

— 

.200 

12 

0.06922 

— 

— 

0.06630 

— 

.200 

13 

0.07476 

— 

— 

0.07137 

— 

.200 

14 

0.08028 

— 

— 

0.07638 

— 

.200 

15 

0.08576 

— 

— 

0.08132 

— 

.200 

16 

0.09121 

— 

— 

0.08621 

— 

.200 

17 

0.09662 

— 

— 

0.09105 

— 

.200 

18 

0.10201 

— 

— 

0.09583 

— 

.200 

19 

0.10736 

— 

— 

0.10056 

— 

.200 

20 

0.11268 

— 

— 

0.10524 

— 

" - "  means  value  is  higher  than  that  given  by  conservative  assumption  of 

independence. 
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Table  Id  Continued . 


2=2.75 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

Trivariat 

P, .  -  P.  . 

U  iJ 

Equicorrelated 
Case  P,  =  P 

P.  .  =  (P 
ij  > 

ij.jj’eT’ 

kl 

ij.jj’eT1 

all  k,2. 

.300 

3 

0.01777 

0.01755 

0.01740 

0.01740 

0.01751 

.300 

4 

0.02363 

0.02334 

0.02304 

0.02295 

0.02326 

.300 

5 

0.02944 

0.02913 

0.02868 

0.02831 

0.02901 

.300 

6 

0.03523 

0.03493 

0.03432 

0.03354 

0.03477 

.300 

7 

0.04098 

0.04072 

0.03996 

0.03867 

0.04052 

.300 

8 

0.04669 

0.04652 

0.04561 

0.04368 

0.04628 

.300 

9 

0.05237 

0.05231 

0.05125 

0.04860 

0.05203 

.300 

10 

0.05802 

— 

0.05689 

0.05342 

0.05778 

.300 

11 

0.06364 

— 

0.06253 

0.05815 

0.06354 

.300 

12 

0.06922 

— 

0.06818 

0.06279 

— 

.300 

13 

0.07476 

— 

0.07382 

0.06735 

— 

.300 

14 

0.08028 

— 

0.07946 

0.07183 

— 

.300 

15 

0.08576 

— 

0.08510 

0.07623 

— 

.300 

16 

0.09121 

— 

0.09075 

0.08056 

— 

.300 

17 

0.09662 

— 

0.09639 

0.08482 

— 

.300 

18 

0.10201 

— 

— 

0.08902 

— 

.300 

19 

0.10736 

— 

— 

0.09315 

— 

.300 

20 

0.11268 

— 

0.09721 

—  — 

.400 

3 

0.01777 

0.01730 

0.01704 

0.01704 

0.01724 

.400 

4 

0.02363 

0.02296 

0.02246 

0.02227 

0.02286 

.400 

5 

0.02944 

0.02863 

0.02788 

0.02726 

0.02847 

.400 

6 

0.03523 

0.03430 

0.03329 

0.03208 

0.03408 

.400 

7 

0.04098 

0.03997 

0.03871 

0.03675 

0.03970 

.400 

8 

0.04669 

0.04563 

0.04413 

0.04127 

0.04531 

.400 

9 

0.05237 

0.05130 

0.04955 

0.04567 

0.05093 

.400 

10 

0.05802 

0.05697 

0.05496 

0.04993 

0.05654 

.400 

11 

0.06364 

0.06264 

0.06038 

0.05409 

0.06215 

.400 

12 

0.06922 

0.06831 

0.06580 

0.05813 

0.06777 

.400 

13 

0.07476 

0.07397 

0.07121 

0.06208 

0.07338 

.400 

14 

0.08028 

0.07964 

0.07663 

0.06594 

0.07900 

.400 

15 

0.08576 

0.08531 

0.08205 

0.06971 

0.08461 

.400 

16 

0.09121 

0.09098 

0.08746 

0.07340 

0.09022 

,400 

17 

0.09662 

— 

0.09288 

0.07701 

0.09584 

.400 

18 

0.10201 

— 

0.09830 

0.08054 

0.10145 

.400 

19 

0.10736 

— 

0. 10371 

0.08401 

0.10707 

.400 

20 

0.11268 

— 

0.10913 

0.08740 

0.11268 

" — "  Means  value  is  higher  than  that  given  by  conservative  assumption  of 
independence . 
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Table  Id  Continued. 


z=2. 75 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

P.  .  ,  =  P. . 
ij  ij 

ijJj’cT’ 

Exact  for 
Equicorrelated 
Case  Pk£  e  P 

all  k.S. 

Trivariate 
P.,,  =  (P..) 

lj  ij 

ij  Jj’eT' 

.500 

3 

0.01777 

0.01691 

0.01653 

0.01653 

0.01683 

.500 

4 

0.02363 

0.02238 

0.02162 

0.02131 

0.02222 

.500 

5 

0.02944 

0.02785 

0.02671 

0.02583 

0.02762 

.500 

6 

0.03523 

0.03333 

0.03180 

0.03013 

0.03302 

.500 

7 

0.04098 

0.03880 

0.03690 

0.03423 

0.03841 

.500 

8 

0.04669 

0.04427 

0.04199 

0.03815 

0.04381 

.500 

9 

0.05237 

0.04974 

0.04708 

0.04193 

0.04920 

.500 

10 

0.05802 

0.05522 

0.05218 

0.04556 

0.05460 

.500 

11 

0.06364 

0.06069 

0.05727 

0.04907 

0.05999 

.500 

12 

0.06922 

0.06616 

0.06236 

0.05246 

0.06539 

.500 

13 

0.07476 

0.07164 

0.06745 

0.05575 

0.07079 

.500 

14 

0.08028 

0.07711 

0.07255 

0.05893 

0.07618 

.500 

15 

0.08576 

0.08258 

0.07764 

0.06203 

0.08158 

.500 

16 

0.09121 

0.08806 

0.08273 

0.06504 

0.08697 

.500 

17 

0.09662 

0.09353 

0.08783 

0.06798 

0.09237 

.500 

18 

0.10201 

0.09900 

0.09292 

0.07084 

0.09777 

.500 

19 

0.10736 

0.10448 

0.09801 

0.07363 

0.10316 

.500 

20 

0.11268 

0.10995 

0.10310 

0.07635 

0.10856 

.600 

3 

0.01777 

0.01632 

0.01579 

0.01579 

0.01621 

.600 

4 

0.02363 

0.02150 

0.02043 

0.02003 

0.02127 

.600 

5 

0.02944 

0.02669 

0.02508 

0.02395 

0.02634 

.600 

6 

0.03523 

0.03187 

0.02972 

0.02761 

0.03140 

.600 

7 

0.04098 

0.03705 

0.03437 

0.03106 

0.03647 

.600 

8 

0.04669 

0.04223 

0.03901 

0.03431 

0.04153 

.600 

9 

0.05237 

0.04741 

0.04366 

0.03740 

0.04660 

.600 

10 

0.05802 

0.05259 

0.04830 

0.04035 

0.05166 

.600 

11 

0.06364 

0.05778 

0.05295 

0.04317 

0.05673 

.600 

12 

0.06922 

0.06296 

0.05759 

0.04587 

0.06179 

.600 

13 

0.07476 

0.06814 

0.06224 

0.04848 

0.06686 

.600 

14 

0.08028 

0.07332 

0.06688 

0.05098 

0.07193 

.600 

15 

0.08576 

0.07850 

0.07153 

0.05341 

0.07699 

.600 

16 

0.09121 

0.08368 

0.07617 

0.05575 

0.08206 

.600 

17 

0.09662 

0.08887 

0.08082 

0.05802 

0.08712 

.600 

18 

0.10201 

0.09405 

0.08546 

0.06022 

0.09219 

.600 

19 

0.10736 

0.09923 

0.09011 

0.06236 

0.09725 

.600 

20 

0.11268 

0.10441 

0.09475 

0.06444 

0.10232 
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Table  Id  Continued 


z=2. 75 


RHO 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

P.  .  =  P. . 
ij  il 

ij.jj'eT' 

Equicorrelated 
Case  Pkl  ,  P 

all  k, 2. 

.700 

3 

0.01777 

0.01546 

0.01476 

0.01476 

.700 

4 

0.02363 

0.02021 

0.01880 

0.01832 

.700 

5 

0.02944 

0.02496 

0.02285 

0.02154 

.700 

6 

0.03523 

0.02971 

0.02690 

0.02447 

.700 

7 

0.04098 

0.03446 

0.03094 

0.02718 

.700 

8 

0.04669 

0.03921 

0.03499 

0.02971 

.700 

9 

0.05237 

0.04396 

0.03904 

0.03208 

.700 

10 

0.05802 

0.04872 

0.04308 

0.03431 

.700 

11 

0.06364 

0.05347 

0.04713 

0.03643 

.700 

12 

0.06922 

0.05822 

0.05118 

0.03845 

.700 

13 

0.07476 

0.06297 

0.05522 

0.04036 

.700 

14 

0.08028 

0.06772 

0.05927 

0.04221 

.700 

15 

0.08576 

0.07247 

0.06332 

0.04397 

.700 

16 

0.09121 

0.07722 

0.06736 

0.04567 

.700 

17 

0.09662 

0.08197 

0.07141 

0.04730 

.700 

18 

0.10201 

0.08672 

0.07546 

0.04888 

.700 

19 

0.10736 

0.09145 

0.07950 

0.05041 

.700 

20 

0.11268 

0.09622 

0.08355 

0.05189 

.800 

3 

0.01777 

0.01417 

0.01332 

0.01332 

.800 

4 

0.02363 

0.01828 

0.01657 

0.01606 

.800 

5 

0.02944 

0.02238 

0.01983 

0.01844 

.800 

6 

0.03523 

0.02649 

0.02309 

0.02057 

.800 

7 

0.04098 

0.03059 

0.02634 

0.02249 

.800 

8 

0.04669 

0.03470 

0.02960 

0.02426 

.800 

9 

0.05237 

0.03880 

0.03285 

0.02588 

.800 

10 

0.05802 

0.04291 

0.03611 

0.02740 

.800 

11 

0.06364 

0.04701 

0.03936 

0.02882 

.800 

12 

0.06922 

0.05112 

0.04262 

0.03076 

.800 

13 

0.07476 

0.05522 

0.04587 

0.03143 

.800 

14 

0.08028 

0.05933 

0.04913 

0.03263 

.800 

15 

0.08576 

0.06343 

0.05238 

0.03378 

.800 

16 

0.09121 

0.06754 

0.05564 

0.03487 

.800 

17 

0.09662 

0.07164 

0.05889 

0.03592 

.800 

18 

0.10201 

0.07575 

0.06215 

0.03692 

.800 

19 

0.10736 

0.07986 

0.06540 

0.03789 

.800 

20 

0.11268 

0.08396 

0.06866 

0.03882 

Trivariate 

P. .  =  (P. . 
iJ  ij 
ij.jj'eT’ 


0.01529 

0.01987 

0.02444 

0.02902 

0.03360 

0.03818 

0.04275 

0.04733 

0.05191 

0.05649 

0.06106 

0.06564 

0.07022 

0.07480 

0.07937 

0.08395 

0.08853 

0.09311 

0.01392 

0.01778 

0.02164 

0.02550 

0.02935 

0.03321 

0.03707 

0.04093 

0.04478 

0.04864 

0.05250 

0.05636 

0.06022 

0.06407 

0.06793 

0.07179 

0.07565 

0.07950 
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Table  Id  Continued. 


z-2.75 


RH0 

Number 

Independence 

Bivariate 

Trivariate 

Exact  for 

P  =  P 

ij  ij 

Equicorre 
Case  P, 

ij.jj'eT* 

kS. 

all  k,  l 

.900 

3 

0.01777 

0.01210 

0.01121 

0.01121 

.900 

4 

0.02363 

0.01517 

0.01339 

0.01293 

.900 

5 

0.02944 

0.01824 

0.01557 

0.01436 

.900 

6 

0.03523 

0.02131 

0.01775 

0.01559 

.900 

7 

0.04098 

0.02438 

0.01993 

0.01667 

.900 

8 

0.04669 

0.02745 

0.02211 

0.01764 

.900 

9 

0.05237 

0.03052 

0.02429 

0.01852 

.900 

10 

0.05802 

0.03359 

0.02648 

0.01933 

.900 

11 

0.06364 

0.03666 

0.02866 

0.02007 

.900 

12 

0.06922 

0.03973 

0.03084 

0.02076 

.900 

13 

0.07476 

0.04281 

0.03302 

0.02141 

.900 

14 

0.08028 

0.04588 

0.03520 

0.02202 

.900 

15 

0.08576 

0.04895 

0.03738 

0.02259 

.900 

16 

0.09121 

0.05202 

0.03956 

0.02314 

.900 

17 

0.09662 

0.05509 

0.04174 

0.02365 

.900 

18 

0.10201 

0.05816 

0.04392 

0.02414 

.900 

19 

0.10736 

0.06123 

0.04610 

0.02462 

.900 

20 

0.11268 

0.06430 

0.04828 

0.02507 

.990 

3 

0.01777 

0.00800 

0.00756 

0.00756 

.990 

4 

0.02363 

0.00902 

0.00814 

0.00796 

.990 

5 

0.02944 

0.01004 

0.00872 

0.00827 

.990 

6 

0.03523 

0.01106 

0.00930 

0.00852 

.990 

7 

0.04098 

0.01208 

0.00987 

0.00872 

.990 

8 

0.04669 

0.01310 

0.01045 

0.00890 

.990 

9 

0.05237 

0.01412 

0.01103 

0.00906 

.990 

10 

0.05802 

0.01514 

0.01161 

0.00920 

.990 

11 

0.06364 

0.01616 

0.01219 

0.00932 

.990 

12 

0.06922 

0.01718 

0.01277 

0.00943 

.990 

13 

0.07476 

0.01821 

0.01335 

0.00954 

.990 

14 

0.08028 

0.01923 

0.01393 

0.00964 

.990 

15 

0.08576 

0.02025 

0.01450 

0.00973 

.990 

16 

0.09121 

0.02127 

0.01508 

0.00981 

.990 

17 

0.09662 

0.02229 

0.01566 

0.00989 

.990 

18 

0.10201 

0.02331 

0.01624 

0.00996 

.990 

19 

0.10736 

0.02433 

0.01682 

0.01003 

.990 

20 

0.11268 

0.02535 

0.01740 

0.01010 

Trivariate  . 

P,  .  =  (P.  .)' 
ij  ij 

ij.jj’eT’ 


0.01183 

0.01452 

0.01732 

0.02005 

0.02284 

0.02556 

0.02835 

0.03104 

0.03383 

0.03652 

0.03932 

0.04206 

0.04474 

0.04755 

0.05024 

0.05302 

0.05571 

0.05850 


means  results  were  unobtainable  due  to  inability  to  integrate  accurately. 
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Table  2A 


Bivariate  Method  Upper  bound  for  Zq  so  that  Pr{^u^X^  i 

where  X  ^  N(0,Z)  and  there  is  a  tree  T  connecting 

such  that  |P  |  s  P  for  all  e  cT. 

x-x-  ij 


(-*oWi>> 

the  X . ,  i=l , . 
1  ’ 


<  .05 


n 

Independence 

P=.l 

P= .  2 

II 

II 

P= .  5 

P= .  6 

4 

2.494 

2.493 

2.491 

2.486 

2.478 

2.467 

2.449 

5 

2.572 

2.572 

2.570 

2.565 

2.558 

2.546 

2.529 

6 

2.635 

2.635 

2.633 

2.628 

2.621 

2.610 

2.593 

7 

2.686 

— 

2.685 

2.681 

2.674 

2.663 

2.646 

8 

2.731 

— 

2.730 

2.726 

2.719 

2.709 

2.692 

9 

2.769 

— 

2.769 

2.765 

2.759 

2.748 

2.732 

10 

2.803 

— 

2.803 

2.800 

2.793 

2.783 

2.767 

11 

2.834 

— 

2.834 

2.831 

2.825 

2.815 

2.799 

12 

2.862 

— 

2.862 

2.859 

2.853 

2.843 

2.828 

13 

2.887 

— 

2.887 

2.884 

2.876 

2.869 

2.854 

14 

2.910 

— 

— 

2.908 

2.902 

2.893 

2.878 

15 

2.932 

— 

— 

2.929 

2.924 

2.915 

2.900 

16 

2.952 

— 

— 

2.950 

2.944 

2.936 

2.921 

17 

2.970 

— 

— 

2.968 

2.963 

2.955 

2.940 

18 

2.988 

— 

— 

2.986 

2.981 

2.973 

2.958 

19 

3.004 

— 

— 

3.003 

2.998 

2.990 

2.976 

20 

3.020 

— 

— 

3.018 

3.014 

3.005 

2.992 

1 - "  means  value  is  higher  than  that  given  by  conservative  assumption 

of  independence. 


rrjr. 


iw 


.yc* 


Table  2B 


it 

Bivariate  Method  Upper  Bound  for  Zq  so  that  Pr{  i  (- 

where  X  ^  N(0,I)  and  there  is  a  tree  T  connecting  the 


such  that 


x  .x . 
i  2 


=  P  for  all  e . .  eT. 


zO0i’zO°i)}  i 
X^  >  i~l » • • • 9 n 


n 

Independence 

P=.l 

P=,2 

P=.  3 

II 

P=  .5 

P-.6 

4 

2.234 

2.233 

2.229 

2.222 

2.211 

2.196 

2.174 

5 

2.319 

2.319 

2.315 

2.309 

2.298 

2.283 

2.261 

6 

2.386 

— 

2.384 

2.378 

2.368 

2.353 

2.331 

7 

2.442 

— 

2.441 

2.435 

2.426 

2.411 

2.390 

8 

2.490 

— 

2.490 

2.484 

2.475 

2.461 

2.440 

9 

2.531 

— 

— 

2.526 

2.517 

2.504 

2.483 

10 

2.568 

— 

— 

2.564 

2.555 

2.542 

2.522 

11 

2.601 

— 

— 

2.597 

2.589 

2.576 

2.556 

12 

2.630 

— 

— 

2.627 

2.619 

2.607 

2.587 

13 

2.657 

— 

— 

2.655 

2.647 

2.635 

2.616 

14 

2.682 

— 

— 

2.680 

2.673 

2.661 

2.642 

15 

2.705 

— 

— 

2.704 

2.696 

2.685 

2.666 

16 

2.726 

— 

— 

2.725 

2.718 

2.707 

2.688 

17 

2.746 

— 

— 

2.746 

2.739 

2.727 

2.709 

18 

2.765 

— 

— 

2.765 

2.758 

2.747 

2.729 

19 

2.783 

— 

— 

2.783 

2.776 

2.765 

2.747 

20 

2.799 

— 

— 

2.799 

2.793 

2.782 

2.765 

.10 


ii 


means  value  is  higher  than  that  give  by  conservative  assumption  of 
independence . 


Table  2B  Continued. 


P= .  7 

P=.  8 

II 

P= .  95 

P= .  99 

4 

2.142 

2.093 

2.008 

1.930 

1.793 

5 

2.228 

2.177 

2.085 

1.997 

1.835 

6 

2.298 

2.246 

2.150 

2.055 

1.873 

7 

2.357 

2.304 

2.206 

2.106 

1.908 

8 

2.407 

2.355 

2.254 

2.151 

1.941 

9 

2.451 

2.399 

2.298 

2.192 

1.972 

10 

2.490 

2.438 

2.336 

2.229 

2.001 

11 

2.525 

2.473 

2.371 

2.263 

2.028 

12 

2.557 

2.505 

2.403 

2.294 

2.054 

13 

2.585 

2.534 

2.433 

2.323 

2.078 

14 

2.612 

2.561 

2.460 

2.350 

2.100 

15 

2.636 

2.586 

2.485 

2.375 

2.122 

16 

2.659 

2.609 

2.508 

2.398 

2.143 

17 

2.680 

2.631 

2.530 

2.420 

2.162 

18 

2.700 

2.651 

2.551 

2.440 

2.181 

19 

2.719 

2.670 

2.570 

2.460 

2.200 

20 

2.737 

2.688 

2.589 

2.478 

2.216 

Table  2C 


Bivariate  Method  Upper  Bound  for  zQ  so  that  X^  i  (~zQ°i ,z0°i) "  — 

where  X  ^  N(0,l)  and  there  is  a  tree  T  connecting  the  X^  i=l,...,n 

such  that  lpx  x  I  =  p  for  a11  eijcT* 

i  j 


n 

Independence 

P=.l 

P= .  2 

P= .  3 

P= .  4 

P-.5 

P= .  6 

4 

1.943 

1.941 

1.935 

1.926 

1.911 

1.891 

1.863 

5 

2.036 

— 

2.033 

2.023 

2.009 

1.989 

1.960 

6 

2.111 

— 

2.109 

2.100 

2.087 

2.067 

2.039 

7 

2.172 

— 

— 

2.164 

2.151 

2.132 

2.104 

8 

2.224 

— 

— 

2.218 

2.205 

2.187 

2.160 

9 

2.269 

— 

— 

2.265 

2.253 

2.235 

2.204 

10 

2.309 

— 

— 

2.306 

2.294 

2.277 

2.251 

11 

2.344 

— 

— 

2.343 

2.332 

2.315 

2.289 

12 

2.376 

— 

— 

2.376 

2.365 

2.348 

2.324 

13 

2.406 

— 

— 

2.406 

2.395 

2.379 

2.355 

14 

2.432 

— 

— 

— 

2.423 

2.408 

2.384 

15 

2.457 

— 

— 

— 

2.449 

2.434 

2.410 

16 

2.480 

— 

— 

— 

2.473 

2.458 

2.435 

17 

2.502 

— 

— 

— 

2.495 

2.480 

2.458 

18 

2.522 

— 

— 

— 

2.516 

2.501 

2.479 

19 

2.541 

— 

— 

— 

2.536 

2.521 

2.499 

20 

2.559 

— 

— 

— 

2.554 

2.540 

2.518 

*» _ "  Means  value  is  higher  than  that  given  by  conservative  assumption  of 

independence . 
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Table  3 


Selected  Trivariate  and  Bivariate  Upper  Bounds  for  so  that 
n 

Prl  u  X . i (-z_o . ,z  z . ) (  <  .10  where  X  N(0,I)  and  there  is  a  tree  T 
x=l_  i  u  l  u  i  — 

connecting  the  X.:  i=l,...,n  such  that  |P  |  =  P  for  all  e  cT. 

l  xixj  ij 


p. . 
13 

Independence 

Bivariate 

Trivariate  (P..,=P..) 

13  13 

e cT ' 

13  33 

Trivariate 

6ij  *33  ' 

.6 

2.32 

2.26 

2.23 

2.25 

.7 

2.32 

2.23 

2.17 

2.22 

.8 

2.32 

2.18 

2.11 

2.17 

.9 

2.32 

2.09 

2.02 

2.07 

ij  13' 


Independence  Bivariate  Trivariate  (P..,*P..)  Trivariate  (P.. ) 

13  13  13  13 


e..  . . , cT ' 

13  ,33 

e.  .  .  .  , 
13  ,33 

2.57 

2.52 

2.49 

2.50 

2.57 

2.49 

2.43 

2.48 

2.57 

2.44 

2.35 

2.42 

2.57 

2.33 

2.24 

2.30 

P..  Independence  Bivariate  Trivariate  (P. .,=P..)  Trivariate  (P . . (P . . ) z) 
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